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EXECUTIVE  SUMMARY 

We  report  on  the  technical  progress  in  the  second  year  of  an  interdisciplinary  experimental,  numer¬ 
ical  and  theoretical  program  to  extend  the  state  of  knowledge  and  understanding  of  the  effects  of 
chemical  reactions  in  hypervelocity  flows.  The  program  combines  unique  experimental  and  compu¬ 
tational  facilities  at  Caltech  to  address  the  key  problems  in  aerothermochemistry  that  arise  from 
interactions  between  the  three  strongly  nonlinear  effects:  Compressibility;  vorticity;  and  chemistry. 
Fig.  1  presents  an  overview  of  the  program. 
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Fig.  1.  Sub-projects,  principal  investigators,  and  main  interactions. 


The  program  has  led  to  the  following  important  new  results  (partial  list): 


•  Clear  understanding  of  the  two  important  parameters  that  define  hypervelocity  flow  over 
spheres.  Closed-form  solution  for  stand-off  distance. 

•  Completion  of  computation  of  hollow-core  compressible  vortex  streets. 

•  First  high-resolution  interferograms  of  shock-shock  interaction  in  hypervelocity  flow.  Detailed 
experimental  and  theoretical  study  shows  that  high-enthalpy  real-gas  effects  do  not  further 
increase  heat  flux  in  type  IV  interaction.  Interferograms  establish  flow  quality  in  T5. 

•  Upgrade  of  S^L  facihty  to  higher  Mach  number.  New  results  on  supersonic  shear  layers  and 
shear- layer /shock- wave  interaction. 

•  Computational  discovery  of  a  flow  field  that  is  sensitive  to  vibration-dissociation  coupling. 

•  Extension  of  the  equilibrium  flux  method  to  a  more  robust,  less  dissipative  form,  and  tests. 

•  Application  of  Schwinger  multi-channel  technique  to  computation  of  collision  cross  section  to 
electrons  of  low  excited  states  of  OH,  NO,  and  002- 

•  Establishment  of  LITA  as  an  accurate  diagnostic  for  gas  properties  over  large  pressure  ranges,  iv 


The  research  has  again  benefited  extensively  from  close  interaction  within  the  program,  including 
regular  research  conferences  and  seminars,  and  through  extensive  communication  with  outside 
groups  through  presentations  at  conferences  and  seminars,  as  well  as  through  consultancies  and  .a 
industrial  participation,  see  individual  projects.  The  progress  made  corresponds  to  more  than  the  — 

half-way  point  according  to  our  plans.  The  research  has  led  to  3  Ph.  D.,  1  Ae.  E.  and  2  M.  S.  _ 

degrees  and  produced  8  new  archival  pubUcations  (-1-7  submitted)  and  8  conference  papers. 
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INTRODUCTION 


Background 

This  section  and  the  following  Statement  of  Work  are  identical  to  the  corresponding  sections 
in  the  first  Annual  Technical  Report  issued  in  May  1994,  see  Candler  et  al,  (1994).  The  repetition 
permits  this  second  report  to  be  read  independently  of  the  first  one.  For  a  complete  picture  of  the 
progress  made  up  to  this  half-way  point  in  the  course  of  this  URI  project,  it  is,  however,  necessary 
to  consult  both  reports. 

At  the  outset  of  this  report  it  is  pertinent  to  quote  three  of  the  opening  paragraphs  of  the 
proposal  that  led  to  the  current  URI  Grant: 

“Transport  to  or  from  space,  aero-assisted  orbital  transfer  operations  and  other  hypervelocity 
flight  maneuvers  necessarily  occur  at  such  high  (near  orbital)  speeds,  that  the  bow-shock  heating  of 
the  components  of  air  causes  them  to  be  dissociated,  to  react  with  each  other,  and,  at  particularly 
high  speeds,  to  be  ionized.  Efficient  air-breathing  propulsion  of  such  high-speed  vehicles  also 
requires  supersonic  combustion  of  the  fuel.  To  a  large  degree,  these  applications  have  motivated 
the  interest  in  the  field  of  study  that  has  been  named  aerothermochemistry  by  von  Karman. 

The  character  of  the  dynamics  of  reacting  gases  is  changed  relative  to  perfect-gas  dynamics  in 
two  important  ways:  First,  the  equilibrium  behavior  is  much  more  complicated,  and  second,  charac¬ 
teristic  times  are  introduced  by  the  fact  that  the  reactions  proceed  at  finite  rates.  These  two  types 
of  highly  nonlinear  gas  imperfections,  interacting  with  the  classical  nonlinear  manifestations  of  gas- 
dynamics,  namely  shock  waves  and  turbulence,  produce  many  new  flow  structures  and  phenomena. 
Some  of  these  are  known  and  qualitatively  understood,  a  few  are  quantitatively  predictable  by 
computation,  aU  are  of  scientific,  and  most  of  practical  importance  to  the  operations  listed  above. 
The  experimental  data  in  this  field  are  very  limited,  because  of  the  difficulties  associated  with 
laboratory  simulation.  This  constitutes  the  most  important  need  in  aerothermochemistry. 

The  proposed  research  work  may  be  grouped  under  the  following  major  headings: 

supersonic  turbulent  mixing  and  combustion 

effects  of  chemical  reactions  in  external  hypervelocity  flows 

reaction  rate  models 

The  philosophy  of  this  proposal  is  to  concentrate  interdisciphnary  effort  on  particular  problems  in 
the  first  two  categories.  The  main  thrust  of  the  proposed  research  is  to  use  the  two  new  facilities  (T5 
hypervelocity  shock  tunnel  and  S^LSupersonic  Shear-Layer  Facility)  for  an  experimental  program 
to  study  and  understand  the  effects  introduced  by  chemical  reactions.  A  parallel  theoretical  and 
computational  effort  will  help  design  the  experiments  and  test  models  and  mechanisms  of  the  gas 
behavior.  The  third  category  pervades  the  whole  program  and  represents  research  in  theoretical 
and  computational  chemistry  and  fluid  dynamics,  as  well  as  comparison  with  experiment.” 
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The  program  was  then  subdivided  into  the  following  four  parts 

I.  Reaction-rate  controlled  shear  flows 

II.  Shock- vorticity  interactions 

III  Supersonic  mixing  and  combustion 

IV  Chemistry  in  non-uniform  flow, 

aU  of  which  involve  complicated  and  non-linear  interaction  between  vorticity,  chemistry  and  com¬ 
pressibility.  By  focusing  the  available  facihties  and  expertise  on  the  main  problems  in  the  field, 
sub-projects  were  determined  as  shown  in  Fig.  1  together  with  the  PFs  and  the  main  interactions 
between  them. 

A  brief  outline  of  the  objectives  was  then  formulated  in  the  following  statement  of  work.  The 
numbering  of  these  sub-projects  corresponds  to  that  of  the  sections  in  this  report  under  which  the 
progress  is  described  in  more  detail. 

Statement  of  Work 

1.1  Obtain  qualitative  information  on  the  structure  of  transitional  and  turbulent  boundary  layers 
at  high  enthalpy.  Determine  the  response  of  a  laminar  boundary  layer  to  disturbances  at  the 
most  strongly  amplified  frequency  for  the  second  mode  instability.  Perform  T5  experiments  on 
a  sharp  slender  cone,  using  differential  interferometry,  laser-induced  fluorescence  and  surface 
instruments. 

1.2  Design  experiments  for  testing  the  limits  of  binary  scaling  by  making  a  numerical  study  of 
blunt  body  flows  in  nitrogen  and  air.  Carry  out  T5  experiments  on  spheres  using  surface 
instrumentation,  differential  interferometry  and  degenerate  four- wave  mixing  (DFWM).  Extend 
the  experiments  to  study  the  high- vorticity  layer  downstream  of  a  bow  shock  associated  with 
the  large  density  rise  caused  by  dissociation. 

1.3  Determine  the  flux  of  reactive  and  inert  species,  and  of  heat  and  momentum  to  the  surface 
of  ablating  bodies  in  hypervelocity  flow.  Carry  out  T5  experiments  using  sampling  mass- 
spectrometric  analysis  and  surface  heat  transfer  and  pressure  measurements.  Formulate  ana¬ 
lytical  models  of  reactive  effects  on  ablation. 

II.  1  Perform  ultra-high  resolution  numerical  experiments  combined  with  nonlinear  analysis  and 
numerics  to  estabhsh  appropriate  initial  conditions  to  elucidate  the  generic  mechanisms  of 
shock- vortex  interactions,  namely,  (1)  interaction  of  a  shock  with  a  vorticity  field  in  the  presence 
of  nonequilibrium  chemistry,  and  (2)  vorticity  generation  in  shock  focusing  with  chemistry. 

II. 2  Determine  the  effects  of  chemical  relaxation  on  the  shock-on-shock  problem.  Demonstrate  the 
magnitude  of  nonequilibrium  effects  on  shock  impingement  heating.  Carry  out  experiments 
with  a  2-dimensional  model  in  T5,  supported  by  local  shock-polar  analysis  and  numerical 
calculations  of  an  ideal  dissociating  gas. 

Ill  Carry  out  S^L  experiments  with  chemically  reacting  and  non-reacting  flows,  using  schlieren, 
spectroscopy,  Rayleigh  scattering,  chemical  effects  and  new  diagnostic  methods,  to  gain  under¬ 
standing  of  the  mechanisms  of  mixing  when  a  shock  interacts  with  a  supersonic  turbulent  shear 
layer,  and  in  the  flow  downstream  of  an  expansion  wave.  Develop  new  related  diagnostics  and 
computational  capabihties. 
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IV.  1  Design  and  carry  out  an  experiment  in  T5  to  provide  data  to  test  chemical  reaction  models  for 

flows  with  vibrational  nonequilibrium.  Use  a  previously  developed  computational  technique 
to  simulate  the  flow  over  proposed  experimental  conflgurations.  Identify  a  configuration  that 
exhibits  significant  effects  of  vibrational  nonequilibrium  on  the  reaction  rates. 

IV. 2  Exploit  computational  methodologies  and  high-performance  and  cost-effective  computing  pro¬ 
vided  by  highly  parallel  and  scalable  supercomputers  to  solve  the  complex  equations  which 
govern  electron- molecule  collision  processes  needed  for  robust  modeling  of  the  chemical  and 
physical  properties  of  hypervelocity  flows. 

IV. 3  Demonstrate  the  effect  of  chemistry  on  the  leeward  shock-vortex  separation  region  of  hyper¬ 
velocity  flow  about  slender  bodies  using  high- resolution  CFD  on  the  Intel  Touchstone-Delta 
computer,  combined  with  dimensional  analysis  and  physical  scaling. 

V. 2  Develop  DFWM  for  measuring  species  concentration  and  temperature  in  hypervelocity  flows. 

concentrate  on  NO-concentrations  because  of  their  importance  in  air  flows.  Buy  the  remaining 
equipment  to  complete  the  set  of  components,  assemble  the  instrumentation,  debug  and  apply 
it  to  T5  flows. 


Current  Status  of  Progress 

In  the  following  paragraphs,  which  mirror  the  statement  of  work,  the  progress  in  the  sub- 

projects  is  summarized. 

1.1  The  experimental  results  of  the  previous  reporting  period  were  further  analyzed,  prepared  and 
submitted  for  archival  publication.  The  experiments  and  equipment  for  sheet  visualization  of 
the  cone  boundary  layer  and  for  disturbance  excitation  at  the  second-mode  frequency  have 
been  prepared.  The  new  experiments  wiU  benefit  from  a  data  acquisition  upgrade  that  was 
implemented  for  an  industrial  user,  giving  much  higher  data  density.  An  existing  chemically- 
reacting  boundary  layer  code  has  been  implemented. 

1.2  Analysis  of  the  results  led  to  a  well  rounded  and  clear  understanding  of  the  two  important 
parameters  that  define  high-enthalpy  flow  over  spheres  in  complex  gas  mixtures  such  as  air 
and  carbon  dioxide,  defining  the  limits  of  binary  scaling  and  giving  a  closed-form  expression 
for  the  shock  stand-off  distance.  In  preparation  for  the  next  phase,  a  conical  nozzle  for  T5  was 
designed,  built  and  partially  calibrated.  This  provides  the  capability  to  operate  satisfactorily 
with  many  different  gases  over  the  whole  enthalpy  range  of  T5.  A  novel  visualization  technique 
similar  to  the  smoke  wire  method  was  demonstrated  and  is  being  further  developed. 

1.3  This  subproject  was  scheduled  to  start  at  the  end  of  the  second  year,  z.e.,  now.  In  the  proposal, 
we  pointed  out  that  it  was  predicated  on  our  ability  to  raise  funds  from  other  sources  to  build 
a  mass  spectrometer.  To  date,  we  have  not  been  able  to  raise  these  funds. 

II.  1  Computation  of  hoUow-core  compressible  vortex  streets  was  completed  and  submitted  for  pub¬ 
lication.  Initial  results  from  extensive  computations  of  a  shock  with  a  contact  discontinuity 
indicate  that  the  early  vorticity  generation  is  not  sgnificantly  influenced  by  dissociative  reac¬ 
tions.  Analytic  work  led  to  an  estimate  for  the  vorticity  generation  by  the  shock/contact-surface 
interaction  and  a  technique  for  computing  similarity  solutions  for  shock/contact  interactions 
with  and  without  chemistry. 
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11. 2  The  completion  of  the  holographic  interferometer  and  the  quality  of  the  heat  flux  gauges  de¬ 
veloped  in  the  previous  reporting  period  came  to  fruition  in  producing  the  first  high-quality 
and  detailed  data  on  shock-shock  interaction  in  hypervelocity  flow.  Results  have  demonstrated 
that  existing  predictions  of  stagnation-point  density  at  high  enthalpy,  owing  to  the  effects  of 
dissociation,  are  not  realized.  Time-resolved  heat  flux  measurements  confirmed  this  observa¬ 
tion  and  agreed  accurately  with  previous  experiments  in  the  low-enthalpy  range.  This  shows 
that  real-gas  effects  do  not  further  increase  the  large  heat  flux  induced  by  shock  impingement. 
A  theoretical  model  was  developed  for  the  influence  of  thermochemistry  on  the  jet  shock  sys¬ 
tem  of  the  type  IV  interaction,  using  the  ideal  dissociating  gas  approximation.  This  leads  to 
the  conclusion  that  peak  heating  arises  from  a  balance  of  the  jet  shock  strengths,  and  that 
this  balance  reduces  the  influence  of  thermochemistry  on  the  flow.  The  new  instrumentation 
developed  for  this  project  enabled  critical  assessment  of  the  flow  quality  in  T5.  High-resolution 
interferograms  of  the  flow  over  blunt  bodies  agree  with  established  theory  within  the  experi¬ 
mental  uncertainty.  The  new  instrumentation  has  also  greatly  benefitted  other  projects  within 
this  URI  program  as  well  as  industrial  users  of  T5. 

Ill  The  (S^L)  combustion  facility  was  further  upgraded  to  permit  operation  at  higher  Mach  num¬ 
ber,  and  experiments  on  supersonic  shear  layers  as  well  as  on  shear-layer/shock- wave  inter¬ 
action  to  be  completed.  Progress  has  been  made  on  the  development  of  a  robust,  unsplit, 
shock-capturing  numerical  simulation  method  for  one-  and  multi-dimensional  chemically  react¬ 
ing  flows  and  detonations.  In  the  Image  Correlation  Velocimetry  project,  the  final  stages  have 
been  reached  of  the  design  and  fabrication  effort  for  the  development  of  new  CCD-imaging 
technology  to  permit  recording  of  image  pairs  with  small  time  separations  (down  to  1  fis). 

IV.  1  The  work  in  the  previous  reporting  period  showed  that  the  flow  over  a  sphere  was  not  suitable 
for  testing  the  vibration-dissociation  coupling  models,  because  it  is  relatively  insensitive  to 
them.  A  computational  search  led  to  the  discovery  of  a  phenomenon  that  is  very  sensitive  to 
the  coupling  model,  and  to  the  design  of  an  experiment  to  test  such  models.  The  configuration 
is  a  double-wedge  model,  in  which  the  shock  on  the  second  wedge  is  detached,  so  that  a 
shock-shock  interaction  results  that  is  very  similar  to  the  Type  IV  interaction  studied  under 
II. 2,  but  with  a  wall-bounded  supersonic  jet.  One  of  the  models  for  the  T5  experiments  has 
been  designed  and  is  being  built.  This  experiment  will  also  benefit  significntly  from  the  data 
acquisition  upgrade. 

IV. 2  Computations  on  the  Intel  Paragon  at  speeds  typically  between  3  and  10  GFLOPS  have  made  it 
possible  to  apply  the  Schwinger  multichannel  method  to  achieve  the  determination  of  electron 
collision  cross  sections  for  three  processes:  1.  Near-threshold  excitation  of  the  lowest  excited 
state  of  OH  (A^S~^).  This  is  important  in  the  dissociation  of  water  vapor.  2.  Impact  excitation 
of  the  A^S"*"  state  of  NO,  which  is  responsible  for  the  7-band  emission  from  the  bow  shock 
of  blunt  body  flows.  3.  Impact  excitation  of  five  electronically  excited  states  of  CO2.  No 
experimental  data  exist  in  any  of  these  cases. 

IV. 3  The  equilibrium  flux  method  was  modified  to  make  use  of  Osher  intermediate  states.  The  new 
scheme  (EFMO)  retains  the  robustness  of  the  EFM  scheme  in  capturing  strong  shocks,  while 
substantially  reducing  its  numerical  diffusivity  for  good  boundary  layer  resolution.  EFMO  was 
incorporated  in  a  fuUy  implicit  version  of  the  viscous  3-D  code  PGP3D.  This  code  was  tested 
extensively  in  flow  over  a  cone  at  zero  and  finite  incidence.  The  ideal  dissociating  gas  model 


-  5- 


has  been  incorporated  into  both  the  implicit  and  the  explicit  versions  of  PGP3D  and  is  being 
tested. 

V.2  This  year,  LITA  has  been  transformed  from  a  promising  technique  into  an  accurate  and  well- 
understood  tool  for  the  study  of  gas  properties.  This  transformation  has  been  accomplished 
through  improvements  in  LITA  signal  modeling  and  a  successful  series  of  LITA  experiments. 
The  experiments  demonstrate  the  accuracy  of  the  improved  model  and  the  ability  to  extract 
accurate  values  of  the  sound  speed  and  thermal  diffusivity  from  LITA  signals.  The  experiments 
also  demonstrate  the  ability  to  obtain  LITA  signals  in  gases  over  a  wide  pressure  range  (~3  Pa- 
14  MPa  at  room  temperature)  and  in  challenging  environments,  e.g.,  CO2  near  its  liquid-vapor 
critical  point,  and  in  flowing,  turbulent  gases. 

A  total  of  7  higher  degrees  were  earned.  8  papers  were  published  or  accepted  for  publication  in 
the  archival  literature,  and  7  have  been  submitted.  9  papers  were  published  as  conference  papers 
or  in  proceedings.  The  work  has  resulted  in  a  large  number  of  presentations  at  various  seminars, 
meetings  and  conferences.  Details  of  these  are  given  in  the  subsections. 

Interactions 

The  interaction  between  the  principal  investigators  and  their  research  groups  within  the  project 
followed  the  diagram  in  Fig.  1.  very  closely.  Much  of  this  interaction  is  forced  by  joint  use  of  facilities 
or  instrumentation,  sharing  of  technical  personnel,  or  through  discussions  that  arise  at  almost 
weekly  interactions  in  the  Fluid  Mechanics  Research  Conference  at  GALCIT.  Interaction  between 
the  separate  Aerothermodynamics  URFs  benefited  from  the  meeting  held  at  Calspan,  Buffalo  in 
1994.  Advantage  was  taken  by  several  PFs  of  attendance  at  conferences  (AIAA,  19thISSW)  and 
meetings  at  AEDC  to  interact  with  other  groups  nationally  and  internationally. 

Significant  industrial  participation  in  the  research  has  begun  through  the  work  of  MetroLaser 
in  the  T5  laboratory.  T5  tests  for  the  European  Space  Agency  and  Japanese  companies  have  also 
been  continued.  A  further  major  international  interraction  is  the  upcoming  20^^  International  Sym¬ 
posium  on  Shock  Waves  organized  to  take  place  at  Caltech  23-27  July  1995.  Some  420  participants 
are  expected  to  attend.  Major  funding  for  this  conference  comes  from  AFOSR. 

Layout 

The  following  chapters  present  the  main  results  of  the  work  in  more  detail,  in  the  order  shown 
above.  Each  subproject  then  lists  the  personnel,  publications,  degrees  earned,  and  interactions 
separately.  A  few  of  the  subprojects  had  received  preliminary  funding  before  the  URI  funding 
began,  so  that  they  are  at  a  more  mature  stage  than  others. 

Reference 

Candler,  G.  V.,  Dimotakis,  P.  E.,  Hornung,  H.  G.,  Leonard,  A.,  Meiron,  D.  L,  McKoy,  B.  V., 
PuUin,  D.  L,  and  Sturtevant,  B.  1994  “Interaction  of  Chemistry,  Turbulence,  and  Shock  Waves  in 
Hypervelocity  Flow,”  AFOSR  Annual  Technical  Report  F49620-93— 1-0338  (GALCIT  Report  FM 
94-2). 
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CHAPTER I 

REACTION-RATE  CONTROLLED  SHEAR  FLOW 


I.l  TRANSITIONAL  AND  TURBULENT  BOUNDARY  LAYERS 


Objectives  and  Status  of  Research 

The  objective  of  this  experimental  investigation  is  to  obtain  qualitative  information  on  the 
structure  of  transitional  and  turbulent  boundary  layers  at  high  enthalpy,  and  to  determine  the 
response  of  the  laminar  boundary  layer  to  disturbances  at  the  most  strongly  amplified  frequency 
for  the  second  mode  instability.  The  experiments  use  surface  heat  transfer  measurements  on  a 
slender,  sharp  cone  and  optical  flow  visuahzation  to  detect  eflFects  of  chemical  reactions  on  the  heat 
transfer,  transition  location,  and  structure  of  the  flow  field. 

The  exploratory  experiments  conducted  during  the  previous  reporting  period  demonstrated 
a  clear  stabilizing  effect  of  high-enthalpy  real-gas  effects  that  increases  with  decreasing  dissocia¬ 
tion  energy,  and  found  that  the  instability  under  the  conditions  tested  occurred  via  the  ToUmien- 
Schlichting  mode.  During  the  last  year,  most  of  the  work  concentrated  on  the  preparation  of  the 
planned  stability  experiment,  the  implementation  of  boundary-layer  software,  and  experiments  at 
higher  enthalpy,  up  to  29  MJ/kg  in  the  laminar  regime  with  nitrogen.  The  gas-dynamical  oscillator 
for  disturbance  generation  proved  to  be  awkward,  and  was  rejected  in  favor  of  the  more  easily 
adjustable  glow  discharge  technique.  The  oscillator  for  this  excitation  is  completed.  The  design  of 
the  tangential-sheet  visualization  in  the  boundary  layer  is  completed,  and  optical  components  are 
being  assembled.  The  extension  of  the  data  acquisition  system  to  more  than  50  test  section  chan¬ 
nels  permits  a  significant  improvement  to  be  made  to  the  data  density  of  the  measurements.  The 
cone  is  being  instrumented  to  accommodate  40  heat  flux  gauges  instead  of  the  14  of  the  previous 
measurements. 

I.  IT  Introduction 

At  sufficiently  high  Mach  number,  theoretical  work  by  Mack  and  others  has  shown  that  in¬ 
stability  to  small  disturbances  in  a  constant-pressure  boundary  layer  occurs  via  a  mode  that  has 
a  characteristic  frequency  far  in  excess  of  the  Tollmien-Schlichting  mode.  Experiments  in  cold 
(perfect-gas)  hypersonic  flows  confirm  this  prediction.  Experiments  in  T5,  obtained  during  the 
previous  reporting  period,  showed  that  the  transition  mechanism  occurs  via  the  T-S  mode.  The 
frequency  of  the  T-S  mode  in  T5  flows  is  typically  150  kHz,  while  the  second- mode  frequency  is  typ¬ 
ically  1  MHz.  The  latter  frequency  is  probably  absent  from  the  noise  spectrum  of  the  facility,  which 
may  be  the  reason  for  the  observed  behavior.  The  experiments  also  showed  that  high-enthalpy  real- 
gas  effects  stabifize  the  boundary  layer,  while  theoretical  results  (Reed  et  a/.,  1992)  indicate  that 
the  second  mode  is  destabilized  by  them.  To  clarify  the  questions  surrounding  these  phenomena,  a 
focussed  stabihty  experiment  is  required. 
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Laminar  constant-pressure  boundary  layers  can  be  accurately  computed  by  boundary  layer 
theory.  The  importance  of  the  chemical  activity  in  the  high-enthalpy  flows  studied  here  must  be 
incorporated  in  the  boundary  layer  calculations.  A  well-established  code  for  this  purpose  is  the 
Boundary  Layer  Integral  Matrix  Procedure  (BLIMP),  which  was  obtained  for  our  use  in  the  URI 
project  from  Air  Force  Wright  Aeronautical  Laboratory.  Comparisons  of  experimental  results  with 
computations  are  required  to  understand  the  influence  of  cheical  reactions  in  the  laminar  boundary 
layer. 


1.1.2  Disturbance  excitation 

High-frequency  oscillations  may  be  generated  in  a  supersonic  flow  by  gas-dynamical  means, 
either  by  a  cavity  in  a  wall,  or  by  a  spiked  blunt  body.  To  obtain  frequencies  in  the  MHz  range,  the 
characteristic  dimension  of  the  blunt  body  or  the  cavity  has  to  be  typically  1  mm  in  a  5  km/s  flow. 
A  spiked  blunt  body  of  this  size  was  tried  in  T5  consisting  of  a  small  blunt  body  (2  mm  diameter) 
with  a  spike.  While  this  method  of  excitation  does  provide  high-frequency  disturbances,  it  is  not 
possible  to  adjust  the  frequency  without  replacing  the  body  or  adjusting  cavity  dimensions. 

The  method  used  by  Kosinovet  a;.(1990)  for  boundary  layer  excitation,  which  employs  a  glow 
discharge  oscillator,  is  more  attractive  from  this  point  of  view.  In  this  technique  a  spark  gap  in  a 
hole  in  the  surface  under  the  boundary  layer  and  almost  flush  with  the  surface  is  excited  from  a 
power-ampUfied  sine  wave  oscillator.  Every  time  the  spark  gap  causes  breakdown,  a  pressure  pulse 
is  emitted.  Kosinov  found  this  to  be  a  suitable  excitation  device.  Our  ten  times  higher  frequency 
requires  the  power  amplifier  to  be  good  up  to  2  MHz.  The  circuit  for  the  excitation  device  has 
been  designed  and  built,  and  is  being  tested. 

1.1.3  Computations  and  comparison  with  experiment 

The  BLIMP  code  solves  the  two-dimensional  or  axisymmetric  boundary  layer  equations  using 
a  fourth-order  Galerkin  scheme  and  a  Newton-Raphson  iteration  procedure  in  a  grid  requiring 
typically  30  points  across  the  boundary  layer.  The  gas  is  modeled  as  a  mixture  of  perfect  gases, 
with  equilibrium,  finite-rate  or  frozen  chemistry.  The  chemical  reactions  are  modeled  with  the 
Arrhenius  equations. 

The  interest  in  using  this  code  is  for  comparisons  with  the  measured  laminar  boundary  layer 
flows.  The  code  also  provides  for  incorporation  of  various  turbulence  models,  for  which  the  tran¬ 
sition  location  is  required  as  an  input.  While  the  comparisons  shown  below  include  computations 
with  various  turbulence  models,  it  is  clear  that  such  models  have  a  much  lower  level  of  scientific 
respectability  than  the  laminar  flow  computations.  Fig.  I.l  shows  results  for  the  case  of  an  inter¬ 
mediate  enthalpy,  intermediate  pressure  flow,  where  the  chemistry  may  be  expected  to  be  frozen. 
It  is  clear  from  the  results  in  the  laminar  range  that  the  agreement  is  satisfactory  there.  The 
discrepancy  between  computation  and  the  similarity  solution  may  be  explained  by  the  difference 
in  evaluating  the  specific  heat  ratio  which  could  be  remedied  by  correcting  the  difference.  The 
differences  between  the  different  turbulence  models  show  that,  even  in  the  perfect-gas  range,  such 
descriptions  of  turbulence  are  inadequate. 

At  higher  specific  reservoir  enthalpy  and  pressure,  where  we  might  expect  equilibrium  behavior 
of  the  gas  in  the  boundary  layer,  the  BLIMP  code  was  applied  in  equilibrium  mode.  The  comparison 
is  shown  in  Fig.  1.2.  In  this  flow  the  boundary  layer  remains  laminar  aU  the  way  to  the  end  of  the 
cone,  because  the  Reynolds  number  is  too  low.  A  significant  difference  is  evident  here,  in  that  the 
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Fig.  LI.  Comparison  of  BLIMP  calculation  with  measurements  of  heat  flux  to  the  slender  cone  in 
the  case  of  a  medium-enthalpy  condition  in  T5,  Open  squares  are  measured  data.  The 
other  symbols  represent  results  for  the  same  conditions  from  BLIMP  computation.  The 
dotted  line  is  the  laminar  boundary  layer  similarity  solution.  The  discrepancy  between  the 
computation  and  similarity  solution  may  be  explained  by  the  difference  in  evaluating  the 
specific  heat  ratio.  Dash-dot  and  dashed  lines  show  application  of  the  turbulence  models  of 
van  Driest  II  and  White  and  Christoph  respectively.  Asterisks,  open  diamonds  and  crosses 
show  BLIMP  computation  with  turbulence  models  of  KendaU,  Bushnell,  and  Cebeci-Smith 
respectively.  In  the  turbulent  layer  computations  the  transition  location  was  taken  from 
the  experimental  results.  The  discrepancy  among  the  turbulence  models  illustrates  the 
inadequacy  of  representing  turbulence  so  simply. 

experiment  shows  a  much  larger  increase  of  the  heat  flux  between  the  frozen  (see  Fig.  LI)  value  and 
the  high-enthalpy  value  than  is  predicted  by  the  computations.  Such  differences  may  be  expected 
to  arise  from  the  fact  that  the  free  stream  in  the  facility  is  frozen  at  a  slightly  higher  dissociation 
fraction  than  the  equilibrium  value,  while  the  calculation  has  an  equilibrium  free  stream  at  the 
same  pressure  and  temperature  as  in  the  facihty.  However,  the  difference  between  the  experimental 
and  computed  heat  flux  incrrease  is  much  larger  than  can  be  explained  by  this  cause.  To  be  able  to 
simulate  the  effect  correctly,  it  is  necessary  to  operate  the  code  in  non-equilibrium  mode  (presently 
being  implemented). 
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Fig.  1.2.  Comparison  of  frozen  and  equilibrium  computations  of  the  laminar-boundary-layer  heat 
flux  with  a  high-enthalpy  experiment.  The  increase  of  the  heat  flux  from  recombination  is 
much  larger  than  predicted. 
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1.2  NONEQUILIBRIUM  AND  VORTICITY  DOWNSTREAM  OF  BOW  SHOCKS 

Objectives  and  Status  of  the  Research 

The  objectives  of  this  sub-project  of  the  program  are  to  design  and  execute  a  series  of  experi¬ 
ments  as  well  as  to  perform  numerical  computations  and  theoretical  analysis,  in  order  to  test  the 
nature  and  limits  of  binary  scaling,  and  to  understand  and  extend  the  reaction-rate  correlation  for 
the  shock-wave  stand-off  distance,  in  reacting  blunt-body  flows.  In  addition,  the  experiments  are  to 
be  extended  to  investigate  the  high-vorticity  layer  downstream  of  the  bow  shock  that  is  associated 
with  the  large  density  rise  caused  by  dissociation.  The  objectives  have  been  extended  relative  to 
those  in  the  original  proposal  because  of  favorable  developments  during  the  reporting  period. 

The  experimental  and  theoretical  work  on  flow  over  spheres  performed  during  the  previous 
reporting  period  with  differential  interferometry  and  heat  flux  measurements  in  T5,  as  well  as  flow 
computations,  were  further  analyzed.  The  result  is  a  well  rounded  and  clear  understanding  of  the 
two  important  parameters  that  define  high-enthalpy  dissociating  flows  over  spheres  in  complex  gas 
mixtures  such  as  occur  in  air  or  carbon  dioxide  flows.  These  lead  to  a  complete  view  of  the  limits 
of  binary  scaling,  and  to  a  closed-form  expression  for  the  shock  stand-off  distance. 

Preparations  for  the  next  stage  of  the  project  were  included  the  design,  manufacture  and  partial 
calibration  of  a  conical  nozzle  for  T5,  and  work  on  the  development  of  a  visualization  technique 
similar  to  the  smoke- wire  method  for  use  in  high-enthalpy  flows. 


1.2.1  Correlation  of  results 

The  results  obtained  in  the  previous  reporting  period  were  further  analyzed  in  order  to  obtain 
a  clearer  picture  of  the  dependence  of  the  flow  field  behind  the  bow  shock  on  a  sphere  on  two  key 
parameters.  The  first  of  these  is  the  reaction  rate  parameter  fl,  the  importance  of  which  emerges 
from  the  following  analysis  of  the  flow  along  the  stagnation  streamline: 

Consider  the  stagnation  streamline  along  the  symmetry  axis  between  the  shock  and  the  stag¬ 
nation  point.  The  momentum  and  energy  equations  for  inviscid  adiabatic  flow  take  the  simple 
forms 

dp  -b  pudu  —  0  =  dh  +  udu 

where  p,  p,  u,  and  h  are  pressure,  density,  velocity  and  specific  enthalpy  respectively.  Thus, 

dp  =  pdh  . 

This  equation  does  not  mean  that  the  entropy  is  constant  along  the  stagnation  streamline,  but 
rather  that  the  only  entropy  change  that  occurs  is  that  associated  with  the  chemical  reaction: 

T ds  —  Ptdci  , 

where  T,  s  are  temperature  and  specific  entropy,  and  the  /ij  and  Cj  are  the  chemical  potentials  and 
mass  fractions  of  the  constituents.  Let  the  caloric  equation  of  state  be  given  in  the  form 


h  =  h{p,p,Ci)  . 
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Since  the  mass  fractions  must  satisfy  the  identity 

n 

Cj  =  1 , 

i=\ 

the  number  of  mass  fractions  that  are  independent  is  one  less  than  the  total  number  n  of  components 
present.  It  is  usually  convenient  to  choose  c\  as  a  dependent  variable  and  the  other  c,  as  independent 
variables.  Thus, 


n  n 

dh  =  hpdp  +  hpdp  +  '^hadci  =  hpdp  +  phpdh  +  y^hadcj 
1=2  1=2 

where  the  subscripts  denote  partial  differentiation.  Solving  for  dp, 


1  —  php  1 

dp  =  — 7 - dh  -  j-^hcidci  . 


t=2 


Note  that  the  coefficient  of  dh  is  related  to  the  frozen  speed  of  sound  ay,  and  that  dh  may  be 
replaced  by  -udu.  Rewriting  the  first  term  on  the  right  of  this  equation  accordingly,  it  becomes 


dp 

P 


v?  du 
aj,  u 


1  -  ^ 
^  ^  2  =  2 


where 

2  _  -hp 
d  hp  ~  1/p  ' 

The  frozen  Mach  number  u/af  after  the  normal  shock  is  typically  0.2  or  less.  This  means  that, 
in  the  absence  of  dissociation,  the  density  is  practically  constant  along  the  stagnation  streamline, 
and,  with  dissociation,  the  density  change  along  the  stagnation  streamline  is  essentially  controlled 
by  the  chemistry: 


(dp)s 


This  approximation  makes  it  possible  to  relate  the  average  density  on  the  stagnation  streamline 
to  the  rate  at  which  energy  is  absorbed  by  the  chemical  reactions  at  the  shock.  This  then  appears 
to  be  the  right  quantity  to  incorporate  in  a  new  reaction  rate  parameter 


For  a  given  gas  mixture. 


O  =  ^  A,  ^ 

”  PsUoo  [hp  dt 


dci 


~  hci  I  /{p. 


\i=2 


oo'^oo)^ 


s 
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which  has  the  physical  significance  of 

~  Energy  absorption  rate  by  chemistry 
Input  rate  of  free  stream  kinetic  energy  * 


The  second  important  parameter  is  the  density  ratio  across  an  equilibrium  normal  shock.  By 
this  parameter,  the  effects  of  the  recombination  rate  are  brought  into  the  analysis  (equilibrium  is 
the  balance  of  forward  and  reverse  reactions.)  This  equilibrium  density  ratio  depends  mainly  on 
the  free  stream  total  enthalpy. 


The  analysis  permits  a  closed  solution  to  be  obtained  for  the  shock  stand-off  distance  in  two 
regimes:  Where  the  equilibrium  density  ratio  does  not  influence  the  result,  and  where,  therefore, 
binary  scaling  works, 


A  = 


h 


-1  + 


(l  + 


1 

2 


(1) 


with  L  —  0.41  and  A  =  Apoofpsd.  At  fi  =  0,  this  goes  to  the  correct  limit  L.  This  is  vabd  when 
the  density  at  the  stagnation  point,  Ph  <  Pe- 


When  Ph  —  Pe^)  binary  scaling  breaks  down,  and  the  appropriate  formula  is 


Again,  this  may  be  seen  to  have  the  correct  limiting  value  Lpsj Pe  when  =  oo. 


(2) 


Equations  1  and  2  are  plotted  in  Fig.  1.3.  The  curves  of  the  two-parameter  family  of  equation 
2,  valid  for  large  fl,  are  nearly  tangent  to  the  single  curve  of  equation  1,  valid  for  small  at  the 
transition  point,  which  is  different  for  different  Ps/Pe* 


The  theory  is  in  good  agreement  with  experimental  results  obtained  over  large  ranges  of  en¬ 
thalpy  in  different  gases,  as  well  as  with  numerically  computed  flows  (see  Wen  and  Hornung,  1995). 


1.2.2  Conical  nozzle 


One  of  the  difficulties  encountered  during  the  experimental  program  was  that  the  contoured 
nozzle  of  T5  is  necessarily  designed  for  one  flow  condition.  The  advantages  of  the  free-piston 
reflected  shock  tunnel  include  the  ease  with  which  the  specific  reservoir  enthalpy  can  be  changed, 
so  that  the  full  range  from  perfect  gas  behavior  to  high-enthalpy  real-gas  behavior  can  be  studied 
in  the  same  facility,  without  changing  the  model  or  test  section  configuration.  However,  in  order 
to  make  fuU  use  of  this  flexibility,  the  contoured  nozzle  has  to  be  used  in  relatively  far  off-design 
conditions  that  lead  to  non-uniform  exit-plane  flow.  This  is  not  the  case  with  a  conical  nozzle, 
since  a  radial  expansion  gives  uniform  flow  on  spherical  surfaces  centered  on  the  source  point  of 
the  flow,  for  any  gas,  whether  reacting  or  perfect.  (At  the  high  densities  and  relatively  small  area 
ratios  in  use  at  T5,  the  nozzle  wall  boundary  layer  is  not  so  important  as  in  nozzles  with  area  ratios 
in  excess  of  200.)  The  penalty  one  pays  by  using  a  conical  nozzle  is  the  associated  xial  gradient 
of  the  flow  properties.  In  some  applications,  such  as  are  to  be  studied  in  this  project  and  in  the 
project  planned  for  the  near  future  by  Candler  (see  section  IV. 1),  the  longitudinal  extent  of  the 
model  is  small,  and  the  axial  gradient  makes  only  little  difference.  Also,  the  conical  free  stream 
can  be  included  in  computations  of  the  flow,  so  that  valid  comparisons  of  measured  and  computed 
flow  fields  can  be  realized. 
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n 


Fig.  1.3.  Plot  of  equations  1  (dashed  line),  and  2  (full  lines)  for  ps/pe  =0.4,  0.5,  0.6,  0.7,  0.8  and 

0.9. 

The  nozzle  was  designed  in  two  parts  so  as  to  be  fully  interchangeable  with  the  existing  con¬ 
toured  nozzle.  The  divegence  angle  of  the  nozzle  was  chosen  at  7  deg.  which  is  a  compromise 
between  sufficient  exit  area  at  a  fixed  nozzle  length  and  small  axial  gradients.  It  was  necessary 
to  make  new  molybdenum  throat  and  nozzle  inserts  for  this  nozzle,  because  the  contoured  nozzle 
starts  with  a  half  angle  of  14  deg. 

A  quick  initial  pitot  calibration  gave  a  flat  distribution  across  the  exit  plane  to  within  the  error 
of  the  pitot-pressure  measurement,  as  expected,  and  in  contrast  to  the  variation  observed  in  the 
off-design  use  of  the  contoured  nozzle  (see  Rousset,  1995). 

1.2.3  Streak-line  visualization 

The  thrust  of  the  project  during  the  next  year  wiU  be  the  investigation  of  flow  fields  in  which 
the  vorticity  generated  by  a  curved  bow  shock  plays  a  role,  this  will  be  particularly  interesting  in 
high-enthalpy  flows,  because  the  vorticity  generated  is  proportional  to  the  density  ratio  across  the 
shock,  and  dissociating  flows  can  give  typically  twice  as  high  a  density  ratio  as  perfect  diatomic 
gas  flows,  and  because  the  shock  curvature  is  also  higher,  which,  for  a  given  body,  further  increases 
the  vorticity,  and  localizes  it  more. 

This  higher  value  and  localization  of  the  vorticity  is  likely  to  cause  the  high-vorticity  layer  to 
become  unstable  and  form  the  kind  of  turbulent  structures  typical  of  a  shear  layer.  To  be  able 
to  see  these,  schlieren  or  interferometric  techniques  will  often  not  be  sensitive  enough.  Also,  a 
streak-fine  visualization  technique  would  be  more  suitable  than  one  that  is  sensitive  to  density  or 
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its  gradient.  For  this  reason  an  attempt  was  made  to  develop  such  a  technique  by  analogy  with  the 
smoke-wire  technique  used  in  subsonic  wind  tunnels. 

Building  on  our  experience  from  the  work  in  the  previous  reporting  period,  an  exploratory 
investigation  was  conducted  in  T5,  in  which  a  wire,  stretched  across  the  nozzle  exit  was  painted 
with  a  weak  aqueous  solution  of  salt.  During  the  flow,  the  salt  is  ablated,  and  dissociates  in  the 
hot  flow.  The  flow  is  then  illuminated  during  the  test  with  parallel  light  tuned  to  the  wavelength 
of  one  of  the  D-lines  of  sodium.  The  sodium  seeding  in  the  flow  absorbs  the  light,  so  that  the  parts 
which  carry  the  sodium  are  opaque,  and  parts  that  are  free  of  sodium  are  transparent. 


Fig.  1.4.  Photograph  of  the  absorption  of  light  by  sodium  seeding  from  a  thin  wire  in  a  high-enthalpy 
flow  in  T5.  Note  the  bending  of  the  wire  under  the  load  of  the  drag  force.  The  flow  is 
opaque  downstream  of  the  wire.  Note  also  the  faint  luminosity  from  the  stagnation  region 
of  the  model. 


The  exploratory  investigation  was  successful  inasmuch  as  it  was  possible  to  demonstrate  that 
the  introduction  of  the  sodium  by  means  of  a  wire  works.  It  also  showed  that  the  sheet  of  sodium 
can  be  made  opaque.  This  may  be  seen  from  Fig.  1.4.  It  is  now  necessary  to  produce  thin  streak 
lines,  and  to  determine  the  degree  to  which  the  wire  affects  the  flow. 
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CHAPTER  II 

SHOCK- VORTICITY  INTERACTION 

11.1  NONEQUILIBRIUM  CHEMISTRY  IN  SHOCK- VORTEX  INTERACTION 

Objectives  and  Status  of  Research 

The  objectives  of  this  research  are  to  elucidate  the  generic  mechanisms  controlling  (1)  inter¬ 
action  of  a  shock  wave  with  a  field  of  vorticity  in  the  presence  of  nonequilibrium  chemistry  and 
(2)  vorticity  generation  mechanisms  in  shock  focusing  with  reaction  chemistry.  Progress  continues 
on  the  computation  of  steady  compressible  vortex  structures;  the  objective  here  is  to  use  these 
solutions  as  special  initial  conditions  for  controlled  numerical  experiments  on  shock  vorticity  in¬ 
teraction.  Work  on  the  computation  of  hollow  core  compressible  vortex  streets  has  now  been 
completed  and  submitted  for  publication.  Work  on  compressible  extensions  of  the  Stuart  vortices 
is  now  also  nearing  completion.  Extensive  computations  of  the  interaction  of  a  shock  with  a  contact 
discontinuity  which  include  dissociation  chemistry  are  also  in  progress.  Initial  studies  show  that 
the  initial  vorticity  generation  due  to  the  shock-contact  interaction  is  not  significantly  affected  by 
the  dissociation  chemistry.  An  analytic  estimate  of  the  deposition  of  vorticity  by  the  shock  has 
been  derived  by  extending  the  shock-polars  analysis  of  Samtaney  and  Zabusky  (1994)  to  include  the 
effects  of  chemistry.  Work  on  the  comparison  of  Whitham  theory  with  numerical  simulation  con¬ 
tinues.  Finally  a  solution  technique  has  been  developed  which  allows  the  computation  of  similarity 
solutions  for  shock-contact  interactions  both  with  and  without  chemistry  effects. 

The  next  stage  of  this  work  wiU  be  the  completion  of  the  comparison  of  the  results  of  Whitham 
theory  with  fuU  numerical  simulation.  In  addition  work  wiU  begin  on  the  use  of  the  steady  vortex 
configurations  computed  above  as  initial  conditions  for  our  flow  simulation  and  the  detailed  study 
of  the  evolving  vorticity  both  in  the  presence  and  absence  of  gas  dissociation. 

11. 1.1  Steady  Compressible  Vortex  Arrays 

As  discussed  in  our  previous  report  part  of  our  study  of  shock-vortex  interactions  progress 
concerns  the  computation  of  steady  compressible  vortex  solutions.  Our  plan  is  to  use  these  solutions 
as  initial  conditions  and  to  then  examine  the  shock- vortex  interaction  resulting  from  the  interaction 
of  an  incipient  shock  with  the  vortex  distribution  both  in  the  presence  and  in  the  absence  of 
chemistry  effects.  Our  work  has  concentrated  on  two  models  of  steady  compressible  vortices.  The 
first  is  an  array  of  hoUow  core  vortices,  an  extension  of  an  approach  initiated  by  Moore  and  PuUin 
(1987).  This  type  of  solution  allows  us  to  use  hodograph  techniques  which  make  it  possible  to  apply 
semi-analytic  techniques  to  effect  the  solution.  This  work  is  now  complete  and  the  results  have 
been  submitted  for  publication.  We  have  shown  that  for  this  particular  compressible  vortex  array 
shock  free  transonic  flow  is  possible  for  certain  choices  of  the  relevant  parameters.  We  have  also 
investigated  the  hmit  of  an  evacuated  vortex  core  and  find  that  we  cannot  obtain  such  solutions 
without  the  formation  of  hmit  fines  indicating  that  shock  free  transonic  flow  may  not  be  possible 
in  this  case. 
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The  second  model  is  obtained  from  a  continuation  of  the  analytical  solutions  of  Stuart  (1967) 
for  a  row  of  steady  incompressible  vortices,  to  finite  Moo-  This  reduces  to  the  solution  of  coupled 
nonhnear  equations  for  the  density  and  stream-functions 


2p^{p^  ^  ”  1)  / dil>\  ( d^\ 

-  1)  ^  ^ 


+  PooP  e' 


-2^/p 


=  0, 


where  p  ~  pj q  —  qjcg^  the  subscripts  “s”  and  “oo”  refer  to  stagnation  conditions  and  conditions 
at  infinity  respectively,  and  ^  is  a  non-dimensional  stream-function.  With  suitable  boundary  con¬ 
ditions  these  equations  can  be  solved  by  finite  difference  techniques  to  yield  steady  ^^compressible 
Stuart  vortex^^  solutions  of  the  Euler  equations.  These  can  then  be  used  as  initial  conditions  for 
the  unsteady  shock-vortex  interaction  codes.  This  work  is  now  well  along.  We  have  again  utilized 
both  Rayleigh-Jantzen  expansions  as  well  numerical  computation  in  order  to  provide  independent 
checks  on  the  results. 


II. 1.2  Numerical  simulation  of  shock  generated  vorticity  with  chemistry 


We  have  extended  our  calculations  of  shock  generated  vorticity  to  hypervelocity  flow  with 
chemistry  effects  included.  We  have  investigated  the  Richtmyer-Meshkov  instability  in  which  a 
shock  impinges  upon  a  contact  discontinuity  for  strong  shocks  i.e.  for  hypervelocity  cases.  To 
model  the  interaction  of  the  flow  with  non-equilibrium  chemical  effects  typical  of  high-enthalphy 
flows,  the  Lighthill-Freeman  ideal  dissociating  gas  model  is  employed.  An  interesting  result  of  these 
investigations  is  that  while  dissociation  causes  significant  changes  in  density  and  temperature,  the 
growth  of  the  perturbations  reduces  neghgibly. 

The  governing  equations  of  motion  are  the  compressible  Euler  equations  given  below. 


Vt  +  Tx{^)  +  Gy{U)  =  S{V), 


where 


V  =  {p,u,v,E,pi,piai,p2a2}  , 

^(U)  =  {pu,pu^  -I-  p,puv,  (E  +  p)u,  piu,  piaiu,  P2a2u}^ , 
^(U)  =  {pv,puv,pv'^  +  p,(E  +  p)v,piv,piaiv,p2a2v}'^, 

rp 

<$(U)  =  {Source  termsp  • 


The  pressure,  p  is  related  to  the  total  energy,  E,  by  p  =  (7  —  1){E  —  ^/>u  •  u). 


The  only  allowable  chemical  reactions  are  the  dissociation  and  recombination  reactions  for  the 
diatomic  species.  Furthermore,  it  is  assumed  that  the  diatomic  species  are  adequately  modeled  by 
the  LighthiU  Ideal  Dissociating  Gas  model.  The  law  of  mass  action  for  the  IDG  model  may  be 
written  as: 


1  -  a; 


^exp  ( 

Pz  \  T 
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where  ttj  is  the  mass  fraction  of  the  dissociated  diatomic  species,  is  the  disscociation  temper¬ 
ature,  pfH  is  the  characteristic  density  and  i  =  1  if  one  dissociating  diatomic  gas  is  present  in  the 
domain  or  i  =  1,2  if  there  are  two  dissociating  gases  in  the  domain.  Note  that  for  the  latter  case, 
the  gases  undergo  dissociation  only  and  do  not  react  with  each  other.  For  non-equilibrium  cases 
the  finite  rate  reactions  proceed  according  to  the  Freeman  model  for  an  IDG  gas  and  are  given 
below: 

^  (<'  -  (^)  -  '’“0  • 

The  function  Cj(T’,  a,)  is  given  by 

C  iT  a  \  =  + 

\Yi 

where  Wi  is  the  molar  mass  (molecular  weight)  of  diatomic  species  i  and  Cu,  C2i,  rjli  and  7?2i 
are  constants  that  can  be  determined  from  experiment.  The  numerical  method  employed  here  is  a 
second  order  accurate  Equilibrium  Flux  Method.  This  technique  now  runs  with  high  efficiency  on 
a  number  of  parallel  architectures  including  the  Intel  Paragon.  We  have  measured  computational 
throughputs  of  roughly  3  Gigaflops  when  this  algorithm  is  run  on  a  100  processor  partition  of  this 
machine. 

It  is  interesting  to  note  that  a  second  order  Godunov  method  seems  to  be  unsuitable  for  strong 
shock  (M  >  10)  cases.  EFM  is  more  diffusive  than  Godunov  type  methods.  However  note  that 
this  feature  of  EFM,  coupled  with  the  fact  that  EFM  does  not  require  the  calculation  of  the  sound 
speed,  is  desirable. 

The  physical  domain  is  a  rectangular  shock  tube;  x  is  the  downstream  direction  and  y  and  2 
are  the  transverse  directions.  The  principal  parameters  are  the  speed  of  the  incident  shock,  Uq,  the 
properties  of  the  gases  across  the  interface  and  the  product  of  the  amplitude  and  wavenumber  of 
the  perturbation,  Ak  =  2‘kA/\.  To  date  we  have  considered  the  following  gas  combinations:  H2- 
N2  (fast-slow)  and  N2-H2  (slow-fast).  The  initial  unshocked  gases  were  maintained  at  a  pressure, 
Po  =  Pb  =  O.latm  and  at  a  temperature  of  Tq  =  Tj  =  298A'.  At  these  temperatures  only  the 
diatomic  species  of  the  gases  were  present. 

Some  representative  results  are  shown  in  the  appended  Figures.  In  Figure  II.l  we  show  density 
contours  for  an  interaction  of  a  Mach  10  shock  with  an  H2-N2  interface  under  frozen  conditions. 
Figure  II.2  shows  the  corresponding  vorticity  for  this  case.  In  Figures  II.3  and  II.4  we  display  the 
density  and  vorticity  respectively  for  the  same  interaction  but  under  equilibrium  conditions. 

II. 1.3  Estimate  of  the  initial  vorticity  using  shock  polars 

It  is  weU-known  that  for  a  perturbation  of  a  fast-slow  interface  where  the  angle  between  the 
incident  shock  and  the  interface  is  constant  and  sufficiently  small  the  Euler  equations  admit  sim¬ 
ilarity  solutions.  In  this  particular  case  the  solutions  correspond  to  constant  states  which  can  be 
obtained  using  the  method  of  shock  polars.  This  technique  has  been  used  by  Samtaney  and  Zabusky 
(1994)to  compute  initial  vorticity  deposition  for  shock  initiated  Richtmyer- Meshkov  instability.  We 
have  extended  these  ideas  to  allow  for  dissociation.  Our  calculations  indicate  that  significant  de¬ 
partures  from  the  frozen  case  occur  at  roughly  Mach  6,  as  indicated  in  Figure  II.5.  However  the 
change  in  the  vortex  sheet  strength  is  small. 
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Figure  II.5.  Local  solution  for  a  equilibrium  dissociating  flow  normalized  by  the  local  solution  for  a 
frozen  case.  The  interface  is  a  sawtooth  with  angle  30°.  The  gases  are  H2-N2.  Here 
C/q,  Co  are  the  shock  speed  and  the  frozen  speed  of  sound  in  the  unshocked  hydrogen,  p2,pt 
are  the  densities  behind  the  reflected  and  transmitted  shock  respectively.  T2,Tt  are  the 
temperatures  behind  the  reflected  and  transmitted  shock  respectively,  and  a  is  the  strength 
of  the  vortex  sheet. 

It  should  be  noted  that  this  initial  estimate  of  the  vorticity  is  only  accurate  at  early  times. 
Subsequent  interactions  due  to  compression  waves  affect  the  circulation  at  the  interface.  In  future 
work  we  will  attempt  to  account  for  these  interactions  so  as  to  obtain  a  more  accurate  estimate  of 
the  growth  rate  of  the  instability  at  later  times. 

II.  1.4  Self-similar  solutions  of  the  compressible  Euler  equations 

We  have  developed  an  iterative  technique  for  computation  of  conically  self-similar  solutions  of 
the  two-dimensional  unsteady  Euler  equations.  The  idea  is  to  search  for  solutions  to  the  Euler 
equations  of  the  form 

U{x,y,t)  = 


where 

^  =  x/t  v  =  y/^ 

Using  this  similarity  transformation  on  the  Euler  equations  leads  to  a  system  of  nonlinear  partial- 
differential  equations 

2U  +  T^  +  Qrj^0 

where 

T  =  {p{u  -  pu{u  -  ()  +  p,  p{u  -  E{u  -  ^)  +  pu}  , 

rri 

Q  =  {p{v  -  ri),  pu{v  -  7?),  pv{v  -7])  +  p,  E{v  -  7?)  -I-  pv}  , 
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which  we  solve  by  fixed  point  iteration.  The  technique  may  be  applied  to  a  fairly  wide  variety 
of  two-dimensional  flows  which  include  the  diffraction  of  a  plane  shock  wave  at  a  corner  and  the 
interaction  of  a  plane  shock  with  a  density  discontinuity  in  the  presence  of  a  wall.  The  similarity 
formulation  includes  flows  with  an  arbitrary  equation  of  state  and  is  therefore  directly  applicable  to 
hypervelocity  flows  with  equilibrium  disscociation/recombination  chemistry  of  the  type  discussed 
above. 

For  the  types  of  initial  conditions  we  have  considered  (such  as  shocks  incident  on  inclined 
contacts)  similarity  solutions  exist  and  this  technique  allows  us  to  compute  them  even  when  the 
subsequent  temporal  evolution  may  exhibit  local  instabilities  of  the  Kelvin-Helmholtz  type  due  to 
the  formation  of  a  contact  discontinuity  (vortex  sheet).  Another  interesting  aspect  of  this  approach 
is  that  much  sharper  resolution  of  the  flow  discontinuities  is  obtained  relative  to  temporal  methods. 
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II.2  SHOCK  WAVE  INTERACTIONS  IN  HYPERVELOCITY  FLOW 

Objectives  and  Status  of  Research 

The  objective  of  this  experimental  study  is  to  determine  the  effects  of  chemical  relaxation  on 
the  shock-on-shock  problem,  and,  in  particular,  to  quantitatively  demonstrate  the  magnitude  of 
nonequilibrium  effects  on  shock  impingement  heating. 

Holographic  interferograms  obtained  for  three  different  test  conditions  and  a  range  of  shock 
impingement  locations  have  demonstrated  that  existing  predictions  of  greatly  increased  stagnation 
density  at  high  enthalpy,  owing  to  the  effects  of  dissociation,  are  not  realized.  Time  resolved  heat 
transfer  measurements  confirmed  this  observation  and  agreed  accurately  with  existing  low  enthalpy 
experimental  investigations.  Thus,  we  have  shown  that  real  gas  effects  do  not  further  increase  the 
large  heat  fluxes  induced  by  shock  impingement. 

The  influence  of  thermochemistry  on  the  jet  shock  system  has  been  modeled  using  the  ideal 
dissociating  gas  approximation.  The  central  conclusion  is  that  peak  heating  arises  from  a  balancing 
of  the  jet  shock  strengths  and  that  this  balance  reduces  the  influence  of  thermochemistry  on  the  flow. 
Real  gas  effects  were  shown  to  be  important  at  lower  Mach  numbers  (<  7.5)  and  for  impinging 
shock  angles  away  from  that  which  produces  peak  heating.  The  model  accurately  reproduces  the 
experimentally  observed  heat  transfer  rates. 

New  instrumentation  developed  for  this  project  has  allowed  critical  assessment  of  the  flow 
quality  produced  by  the  T5  facility.  Measurements  of  the  flow  over  blunt  bodies  without  shock 
impingement  agree  with  established  theories  within  experimental  uncertainty.  This  confirms  the 
validity  of  the  current  data.  This  work  is  reported  in  full  by  Sanderson  [9]. 


II. 2.1  Introduction 

The  interaction  of  a  weak  oblique  shock  with  the  strong  bow  shock  ahead  of  a  blunt  body  in 
supersonic  flow  produces  extreme  heat  transfer  rates  and  surface  pressures.  Although  the  problem 
has  been  studied  extensively  in  low  enthalpy  flows  (Edney  [3]  [4]),  the  influence  of  high  enthalpy 
real  gas  effects  remains  unknown.  Edney  [3]  observed  and  classified  six  distinct  interaction  regimes 
that  are  known  as  types  I— VI.  Severe  heating  occurs  for  the  type  IV  flow  where  the  interaction 
causes  the  formation  of  a  supersonic  jet  that  impinges  on  the  surface  of  the  body.  Local  heat 
transfer  rates  and  surface  pressures,  where  the  type  IV  jet  meets  the  body,  typically  exceed  the 
undisturbed  stagnation  point  values  by  a  factor  of  ten. 

Shock  loci  techniques  developed  for  the  study  of  Mach  reflection  problems  (Courant  k  Friedrichs 
[2])  give  the  solutions  near  the  shock  interaction  points  by  mapping  the  problem  into  the  pressure- 
flow  deflection  angle  {p-6)  plane.  This  method  of  calculation  allowed  Edney  [3]  to  study  the  variation 
of  pressure  at  the  jet  impingement  point  for  a  perfect  gas.  The  results  of  such  a  computation  are 
shown  in  figure  II. 6  for  typical  values  of  the  parameters.  Using  the  stagnation  point  boundary  layer 
solution  of  Cohen  &  Reshotko  [1],  Edney  [3]  correlated  the  heat  transfer  with  the  peak  pressure 
and  body  diameter  to  jet  width  ratio,  Df  Ax,  in  the  form. 


Here  the  peak  heat  transfer  rate  and  pressure  are  normalized  with  respect  to  the  undisturbed 


4) 

^0  /  peak 


PoJ  peak 


(1) 
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Figure  II. 6.  Ratio  of  jet  impingement  pressure  to  undisturbed  stagnation  point  pressure  as  a  function  of 
incident  shock  flow  deflection  angle,  6i,  and  ratio  of  specific  heats,  7  (after  Edney  [3];  Figure  7-9).  Free 
stream  Mach  number  is  —  10. 


stagnation  point  values.  Given  equation  1,  figure  II. 6  indicates  that  reductions  of  the  ratio  of 
specific  heats  at  high  enthalpy  should  produce  greatly  increased  heat  transfer  rates.  The  heat 
transfer  intensification,  iqlqo)peak-'  increases  strongly  with  Mach  number.  Although  the  ratio 
of  specific  heats  is  only  a  crude  measure  of  the  influence  of  equilibrium  thermochemistry,  the 
magnitude  of  the  effect  motivates  the  current  experiments. 


II. 2. 2  Scope  of  the  current  study 

The  objective  of  the  present  work  was  to  determine  the  quantitative  effects  of  thermochemistry  on 
the  shock  impingement  phenomenon.  The  current  work  consisted  primarily  of  a  series  of  experi¬ 
ments  for  both  perfect  and  dissociating  gases  that  utilized  the  capabilities  the  GALCIT  T5  free 
piston  shock  tunnel.  The  flow  topology  was  simphfied  by  studying  the  nominally  two  dimensional 
flow  about  a  cylinder  with  a  coplanar  impinging  shock  wave. 

A  significant  preliminary  objective  of  this  project  was  to  develop  effective  instrumentation 
for  use  in  the  hostile  flow  environment  of  hypervelocity  shock  tunnels.  High  resolution  holographic 
interferometry  was  used  to  investigate  changes  in  the  flow  structure  as  the  location  of  the  impinging 
shock  wave  was  varied.  Fast  response  heat  transfer  gauges  were  developed  to  provide  time  resolved 
measurements  of  the  model  surface  temperature.  The  application  of  these  developments  to  achieve 
the  fluid  mechanical  goals  of  the  project  is  reported  in  the  following  sections. 
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II. 2. 3  Description  of  experimental  apparatus 

The  experimental  arrangement  is  shown  in  figure  11.7.  High  temperature  and  pressure  gas  was 


Figure  II. 7.  Layout  of  apparatus  in  the  test  section  of  T5. 

generated  in  the  nozzle  reservoir  by  the  T5  hypervelocity  shock  tunnel.  Free  stream  conditions 
are  listed  in  table  ILL  The  cylindrical  model  (<?!>40  mm  x  180  mm)  was  positioned  in  the  free  jet 


Test  condition 

A 

B 

C 

Nozzle 

Pressure 

14.6 

22.3 

28.3 

reservoir 

(MPa) 

±0.7 

±1.3 

±1.9 

conditions 

Temperature 

(K) 

3210 

7550 

8960 

Enthalpy 

(MJ/kg) 

3.88 

12.0 

19.1 

Test  section 

Velocity  (m/s) 

2540 

4450 

5350 

conditions 

Density  (kg/m^) 

0.0218 

0.0155 

0.0157 

Pressure  (kPa) 

1.03 

5.48 

11.4 

N  concentration  (kg/mole) 

9.9  X  10"® 

6.9  X  10“^ 

3.65  X  10° 

Moo 

9.9 

6.3 

5.3 

Re  per  mm 

5500 

1600 

1350 

Table  ILL  Summary  of  the  free  stream  conditions.  Results  for  condition  B  are  not  reported  here. 

test  section  through  which  the  test  gas  exhausted  into  an  evacuated  dump  tank.  A  large  shock 
generator  was  partially  inserted  into  the  uniform  nozzle  core  flow  that  extended  upstream  of  the 
nozzle  exit  plane.  Variations  in  the  shock  impingement  flow  fleld  were  investigated  by  vertically 
translating  the  model  with  respect  to  the  exit  of  the  shock  generator.  Optical  access  for  holographic 
interferometry  was  available  through  windows  that  were  mounted  to  flanges  on  the  sides  of  the  test 
section.  Thermocouple  gauges  were  used  to  record  the  surface  temperature  at  24  circumferential 
locations  around  the  model  forebody. 
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II. 2. 4  Experimental  results 

Figures  II. 8  &  II.9  show  holographic  interferograms  and  heat  transfer  distributions  obtained  for 
the  type  IV  interaction  at  low  and  high  enthalpy  respectively.  The  interferograms  exhibit  the 
characteristic  features  of  the  type  IV  flow.  Two  mutually  inverted  A-shocks  are  observed  along 
with  a  terminating  strong  shock  and  a  region  of  elevated  density  behind  it.  Development  of  the 
holographic  interferometer  was  successful  in  resolving  the  internal  structure  of  the  type  IV  jet.  The 
circumferential  locations  of  the  maximum  instantaneous  heat  fluxes,  that  are  indicated  in  figures 
II.8  &  11.9,  coincide  with  the  corresponding  fringe  shift  maxima.  The  peak  heat  fluxes  recorded 
during  the  entire  sequence  of  shots  at  each  test  condition  are  listed  in  table  II. 2. 


Test  condition 

A 

c 

Peak  measured 

Qbody 

Instantaneous 

15.2 

8.0 

100  fis  averaged 

13.5 

6.0 

Predicted 

9body 

Frozen 

14.1 

5.07 

Nonequilibrium 

- 

4.37 

Equilibrium 

- 

7.23 

Table  II.2.  Comparison  of  measured  and  predicted  heat  transfer  intensification  for  type  IV  flows  with 
cylindrical  model. 

The  conclusion  of  the  experimental  investigation  was  that  the  stagnation  density  and  heat 
flux  intensification  decreased  at  high  enthalpy.  Care  must  be  exercised  in  the  interpretion  of  this 
trend  since  the  Mach  number  of  the  flow  in  T5  decreased  as  the  stagnation  enthalpy  was  increased 
(Table  II.l).  The  net  result  was  that  the  influence  of  real  gas  effects  was  observed  to  be  much  less 
than  that  predicted  by  the  variable  7  model  (Figure  II.6). 

In  addition  to  these  measurements  of  the  influence  of  real  gas  effects,  the  high  bandwidth  of  the 
heat  flux  gauges  permitted  examination  of  the  temporal  characteristics  of  the  heat  transfer  signals. 
Sanderson  [9]  discusses  observations  of  coherent  unsteady  fluctuations  of  the  type  IV  jet  that  are 
consistent  with  oscillations  observed  in  recent  numerical  simulations  (Gaitonde  [7],  Zhong  [10]). 


II.2.5  Model  of  the  influence  of  real  gas  effects 

Modeling  of  the  influence  of  real  gas  effects  followed  the  pattern  of  the  existing  perfect  gas  model  of 
the  type  IV  field  (Edney  [3]).  Thermochemical  effects  were  investigated  in  four  conceptual  stages; 

1.  Normal  shock  solutions  were  obtained  using  the  ideal  dissociating  gas  model  due  to  LighthiU 

[8]. 

2.  Extension  to  oblique  shock  waves. 

3.  Mapping  to  the  p-6  plane  to  obtain  solutions  for  the  A-shock  points. 

4.  Modeling  of  the  impingement  heating  using  the  equilibrium  reacting  flow  similarity  solution 
of  Eay  &  Riddell  [5] . 


Figure  IL8.  Shot  T5-850;  condition  A. 

Sanderson  [9]  shows  that  the  heat  transfer  intensification  may  be  expressed  in  the  form; 


Qbody  ^ 


D 


Cl,jtC2jet 


=  func{Poo,  Ho^,  pd,oico,  ^i)- 
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The  parameters  are  Poo  =  PcolPooulo  ~  Hq  is  the  stagnation  enthalpy  normalized  with 

respect  to  the  dissociation  energy  of  the  gas,  pd  is  the  equilibrium  constant,  Ooo  is  the  dissociated 
mass  fraction  of  the  free  stream  and  l3i  is  the  angle  of  the  impinging  shock  wave.  Since  this  is  a 
local  solution  the  geometrical  terms  on  the  left  hand  side,  that  specialize  the  result  for  a  particular 
configuration,  remain  indeterminate;  where  AxfD  is  the  jet  width  to  body  diameter  ratio.  The 
various  similarity  coefficients,  (^,  are  determined  from  the  similarity  solution  of  Fay  &  Riddell  [5] 
and  the  Newtonian  solution  for  the  velocity  gradient  at  the  outer  edge  of  the  boundary  layer. 

The  behavior  of  the  local  solution  2  with  stagnation  enthalpy  is  illustrated  in  figure  II.  10  for 
typical  values  of  the  remaining  parameters.  Three  solutions  were  obtained  in  order  to  investigate 
the  influence  of  nonequilibrium  thermochemistry; 

1.  The  frozen  solution  that  applies  when  the  reaction  rates  are  slow. 

2.  The  equilibrium  solution  that  applies  when  the  reaction  rates  are  fast. 

3.  A  nonequilibrium  approximation  whereby  the  weak  oblique  waves  at  the  A-points  are  frozen 
and  the  strong  bow  shock  and  terminating  jet  shock  reach  equilibrium. 

The  validity  of  the  approximate  nonequilibrium  solution  was  assessed  by  computing  the  reaction 
rates  (Freeman  [6])  immediately  downstream  of  the  translational  discontinuities.  Figure  II.IO  shows 
that  equibbrium  real  gas  effects  increase  the  heat  transfer  intensification  with  respect  to  the  frozen 
solution.  This  occurs  because  of  increased  density  rise  across  the  weak  oblique  jet  shocks.  In  the 


Figure  II.IO.  Influence  of  dimensionless  stagnation  enthalpy,  Ho^,  Mach  number,  ~  1/Pto  and  reaction 

C2  ^ 

rate  on  type  IV  heat  transfer  intensification,  .  . curves  are  the  frozen  solutions.  — 

•  •  —  curves  are  the  equilibrium  solutions. - curves  are  the  nonequilibrium  approximate  solutions.  The 

parameters  are  Poo  =  0.002,  0.010,  0.025,  pd  =:  1  x  10^,  aoo  =  0  and  /?i  =  16®. 
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nonequilibrium  approximation,  where  the  jet  shocks  remain  frozen,  the  heat  transfer  intensification 
decreases  slightly  with  respect  to  the  frozen  solution.  Within  the  scope  of  the  approximation,  non¬ 
monotonic  behavior  with  the  reaction  rate  is  indicated.  Further  analysis  (Sanderson  [9])  demon¬ 
strates  that  peak  heating  conditions  for  a  perfect  gas  minimise  the  influence  of  real  gas  effects 
for  equilibrium  flows  and  promote  nonequilibrium  effects  by  depressing  the  reaction  rate  in  the 
supersonic  jet. 

Specific  predictions  of  the  model  for  the  test  conditions  listed  in  table  II.l  are  compared  with 
the  experimental  data  in  table  IL2.  The  extent  of  chemical  reactions  at  the  low  enthalpy  condi¬ 
tion  A  was  negligble  and  so  only  the  frozen  solution  is  shown.  Reaction  rates  for  condition  C  were 
sufficiently  high  that  departures  from  equilibrium  were  expected  to  be  slight.  The  calculated  values 
accurately  reproduce  the  experimental  heat  transfer  measurements. 
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CHAPTER  III 

SUPERSONIC  SHEAR-FLOW  MIXING  AND  COMBUSTION 


Objectives  and  status  of  the  research 

The  objective  of  this  part  of  the  research  effort  is  to  gain  an  understanding  of  mixing  and 
combustion  processes  in  supersonic  and  compressible  shear  flows,  including  the  effects  of  shock- 
wave  interactions  with  turbulent  mixing  regions.  This  part  of  the  effort  includes  experimental 
investigations  in  supersonic  mixing-layer  flows,  with  supersonic/subsonic  freestreams,  as  weU  as  all- 
hyperbolic  (supersonic/supersonic)  freestreams,  a  numerical  simulation  effort  of  chemically-reacting 
flows  and  detonations,  and  the  development  of  image  correlation  velocity-measurement  techniques 
and  instrumentation  that  would  be  applicable  to  both  subsonic  and  supersonic  flows.* 

Progress  during  the  reporting  period  has  been  made  in  all  three  areas,  with  a  major  upgrade 
of  the  Supersonic  Shear  Layer  (S^L)  Combustion  facility,  to  allow  higher  Mach  number  flows  to 
be  attained,  and  experiments  on  supersonic  shear  layers  and  shear-layer /shock- wave  interactions 
completed;  the  development  of  a  robust,  unsplit,  shock-capturing,  numerical  simulation  method 
for  the  calculation  of  one-  and  multi- dimensional  chemically-reacting  flows  and  detonations;  and 
further  developments  in  Image  Correlation  Velocimetry.  We  are  presently  in  the  final  phases  of  a 
design  and  fabrication  effort  to  develop  new  CCD-imaging  technology  to  permit  the  recording  of  a 
pair  of  images  in  quick  succession  (as  close  as  1  /is  apart).  Several  of  the  new  CCD-imaging  devices 
have  tested  favorably  on  the  wafers  recently  returned  to  us  from  the  silicon  foundry. 

In  the  next  stage  of  this  effort,  direct  measurements  of  mixing  in  all-hyperbohc  shear  layers  wiU 
be  undertaken,  as  weU  as  the  development  of  a  modified  test  section  for  the  S^L  facility  to  study 
wall-bounded  mixing  and  combustion  in  supersonic  flows;  the  numerical  simulations  of  chemically- 
reacting  flows  and  detonations  wiU  be  extended  to  unsteady,  two-dimensional  flows;  and  the  first 
trials  of  the  new  image-pair  recording  CCD’s  wiU  be  conducted  in  a  laboratory  environment. 

II. 2.1  Bi-supersonic  shear-layer  flows 

As  part  of  our  continuing  experimental  investigations  of  the  behavior  of  all-hyperbohc,  fuUy- 
developed,  turbulent,  mixing-layer  flows,  we  have  studied  two  bi-supersonic  flows,  at  low  and 
moderate  compressibility  levels.  An  example  of  the  low-compressibihty  flow  (nominal  nozzle-exit 
conditions:  Mi  =  1.5  [N2],  M2  —  1.2  [N2])  is  shown  in  Fig.  III.l.  Evident  in  this  schheren  flow  visu- 
ahzation  is  a  near-normal  stationary  compression  shock  in  the  bottom  freestream,  located  shghtly 
downstream  of  the  nozzle-exit  plane.  The  (obhque)  Mach  wave  incident  on  the  normal-wave  front 
is  seen  to  terminate  there,  as  demanded  by  the  downstream  (shock-processed)  subsonic  flow.  The 
flow  is  observed  to  be  supersonic  a  short  distance  downstream  of  this  wave.  As  indicated  at  this 
downstream  location,  an  oblique  Mach  wave  is  refracted,  from  the  top  freestream,  through  the  mix¬ 
ing  layer,  into  the  bottom  freestream.  Refractions  of  this  type  are  only  possible  in  aU-hyperboHc 
(as  measured  in  the  Gahlean  reference  frame  of  the  wave)  flow. 

*  The  latter  two  efforts  are  cofunded  under  an  AFOSR  Air-Breathing  Propulsion  program,  Grant  No.  F49620-92- 

J-0290. 


Fig.  III.l.  Schlieren  (oblique  knife  edge)  photograph  of  Mi  =  1.5  [N2],  M2  —  1.2  [N2]  (nominal 
nozzle-exit  conditions)  mixing-layer  flow. 

The  bottom  freestream  can  be  seen  to  be  accelerated  back  to  a  supersonic  state  by  an  effec¬ 
tive  freestream  area  change  produced  by  the  displacement  thicknesses  of  the  mixing  layer  and  the 
bottom-guidewaU  boundary  layer,  in  addition  to  the  weak  curvature  of  the  shear  layer.  This  inter¬ 
action,  dubbed  the  “aerodynamic  Laval  nozzle”,  was  also  described  in  the  discussion  of  a  different 
mixing-layer  flow  in  our  previous  Annual  Technical  Report  (Candler  et  al.  1994).  The  significant 
influence  even  small  changes  in  the  effective  freestream  area  has  on  these  flows  is  noteworthy.  This 
sensitivity  is  largely  traceable  to  the  the  Mach- number /area  relationship,  e.^.,  Fig.  5.2  of  Liepmann 
&  Roshko  (1957),  for  quasi-one-dimensional  flow  in  the  transonic  Mach  number  regime. 

A  second  example  of  low-compressibility  flow  behavior  is  shown  in  Fig.  III.2.  The  non-linear 
computer-control  system  for  the  Supersonic  Shear  Layer  Facility  (c/.  Hall  1991),  permits  extremely 
precise  control  of  the  freestream  stagnation  pressures,  allowing  for  considerable  flexibility  in  setting 
particular  test-section  flow  conditions.  Exploiting  this  ability,  a  flow  configuration  was  designed, 
at  a  slightly  lower  (5%)  top-freestream  stagnation  pressure  than  the  flow  discussed  above.  This 
change  in  the  top-freestream  nozzle-exit  condition  significantly  modifies  the  interaction  between  the 
two  freestream  flows  and  the  fuUy- developed  turbulent  mixing  layer.  Notably,  the  normal  shock 
observed  in  the  previously-discussed  (Fig.  III.l)  flow  is  now  “unfolded”  into  a  Mach  reflection  of 
the  incident  oblique  shock  from  the  bottom  guidewaU.  Visualizations  recorded  in  a  set  of  runs 
of  decreasing  stagnation  pressures  exhibit  a  gradual  shortening  of  the  Mach  stem,  until  regular 
reflection  is  attained. 

These  flow  configurations,  z.e.,  two  moderately-supersonic,  acousticaUy-similar  gas  streams,  pro¬ 
vide  us  with  an  important  diagnostic  capability.  In  these  flows,  turbulence-generated  density  gradi¬ 
ents  are  relatively  weak  in  comparison  to  the  wave-generated  gradients.  As  observed  in  the  schlieren 
photographs,  e.g.^  Fig.  III.l,  this  discrepancy  in  the  two  types  of  gradients,  in  addition  to  the  trans¬ 
mission  of  the  waves  through  the  mixing  layer,  provides  an  instantaneous  visualization  of  the  details 
of  the  wave-turbulence  interactions.  Given  that  many  of  the  these  waves  are  well-approximated 
as  infinitesimal  disturbances,  these  interactions  can  be  viewed  as  passive  measurements  of  the 
instantaneous,  fully-developed  turbulent  Mach  number  field. 

We  have  also  investigated  some  bi-supersonic  mixing-layer  flows  at  moderate-compressibility 
conditions.  A  schlieren  visualization  of  this  M\  =  1.5  [He],  M2  =  1.2  [N2]  (nominal  nozzle-exit 
conditions)  mixing-layer  flow,  is  reproduced  in  Fig.  III. 3.  The  compressibility  of  this  flow,  as 


Fig.  III.2.  Schlieren  (oblique  knife  edge)  photograph  of  Mi  =  1.5  [N2],  M2  =  1.2  [N2]  (nominal 
nozzle  exit  conditions)  mixing-layer  flow.  Visualization  recorded  at  a  slightly  lower  top- 
freestream  nozzle-exit  stagnation  pressure  than  that  of  the  Fig.  IILl  flow. 


measured  by  the  total  convective  Mach  number, 
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is  Me  ^  0.7,  in  contrast  to  the  low-compressibility  flows  (Me  ~  0.1)  discussed  above.  A  consequence 
of  the  increased  compressibility  level  is  the  presence  of  two  different  pressure- wave  systems  in  this 
flow:  a  lab-frame-stationary  system  generated  by  the  imposition  of  the  mixing-layer  displacement 
thickness  when  the  two  freestream  flows  first  interact  (extreme  left-hand  side  of  Fig.  IIL3)  and  a 
traveling-wave  system  in  the  bottom  freestream.  Significantly,  no  such  traveling  wave  system  was 
observed  in  the  low-compressibility  flows. 


As  noted  earlier,  the  stationary-wave  refractions,  through  the  mixing  layer,  observed  in  this  flow 
are  disallowed  if  one  of  the  freestreams  is  elliptic,  in  the  Galilean  frame  of  the  wave.  These  finite- 
strength  wave-turbulence  interactions  do  not  appear  to  have  a  significant  effect  on  the  growth  rate 
of  the  layer,  consistent  with  the  observations  of  Shau  et  al  (1993).  The  insensitivity  of  the  mixing- 
layer  growth  rate  to  these  wave-turbulence  interactions,  in  this  high-density-ratio,  {P2I Pi  ~  5)-case, 
is  interesting,  considering  the  (net)  baroclinic  vorticity  (circulation)  that  such  refractive  interactions 
can  generate. 


The  system  of  traveling  waves,  present  in  the  freestream  flow,  are  of  the  same  character  as 
those  observed  originally  by  Lowson  and  OUerhead  (1968),  Tam  (1971),  and  Oertel  (1979),  in 
supersonic,  turbulent,  jet  flows,  and  more  recently  by  HaU  et  al.  (1993),  in  a  supersonic,  mixing- 
layer  flow.  As  discussed  by  these  authors,  these  waves  appear  to  be  associated  with  large-scale 
turbulent  structures,  convecting  supersonically  with  respect  to  (in  this  case)  the  bottom  (low- 
speed)  freestream.  Defining  convective  Mach  numbers  of  these  features  with  respect  to  each  of  the 
freestreams,  i.e.. 
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where  Uq  is  the  lab-frame  convection  velocity  of  these  features,  we  deduce  Md  ^  0.2  and  Mc2  ^  1-6. 
This  Mast-mode”  convection  velocity,  as  observed  in  a  bi-supersonic  mixing-layer  flow,  is  at  variance 
with  the  convection- velocity  stream-selection  rule  proposed  by  Papamoschou  (1989).  Additional 
exceptions  to  this  proposed  rule  are  found  in  the  experiments  of  Martens  et  al.  (1994)  as  well  as 
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in  the  Papamoschou  data  (1989).**  These  observations  are  at  a  lower  Me,  however,  than  this  case, 
where  the  effect  of  compressibility  appears  to  be  less  significant,  ie.,  Mq\  Mc2  »  Me  for  small 
Mcli  Mc2- 


Fig.  IIL3.  Schlieren  (vertical  knife  edge)  photograph  of  Mi  =  1.5  [He],  M2  =  1.2  [N2]  (nominal 
nozzle-exit  conditions)  mixing-layer  flow. 

This  effort  is  part  of  the  Ph.D.  research  of  Michael  Slessor. 

II. 2. 2  Chemically- reacting,  compressible-flow  simulations 

As  part  of  our  effort  to  accurately  and  efficiently  compute  chemicaUy-reacting  compressible 
flows,  we  have  developed  a  new,  unsplit,  algorithm  for  the  solution  of  the  governing  equations. 
Conventional,  i.e.,  split,  algorithms  employ  separate  numerical  steps  for  the  gas  dynamics,  ignoring 
the  contribution  of  chemical  reactions,  and  treat  the  chemicaUy-reacting  source  terms  in  a  separate 
step,  artificiaUy  decoupUng  the  dynamics.  As  a  consequence,  the  application  of  such  methods  to 
chemicaUy-reacting  flows  maybe  questionable,  depending  on  the  details  of  the  reacting  flow,  because 
of  the  extremely  large  range  of  (coupled)  spatial  and  temporal  scales  encountered  in  such  problems, 
as  weU  as  the  very  large  number  of  unstable  modes  that  characterize  such  systems,  leading  to  a 
multipUcity  of  both  Unear  and  nonUnear  instabiUty  paths.  The  new  algorithm  eUminates  these 
difficulties  by  integrating  the  complete  conservation  equations,  including  the  contributions  by  the 
chemical  reactions,  in  a  single,  fuUy-coupled,  step. 

The  new  algorithm  is  an  extension  and  generaUzation  of  the  scheme  developed  at  Caltech,  under 
AFOSR  sponsorship,  for  non-reacting  flows  (Lappas  1993,  Lappas  et  al.  1994),  finds  surfaces  in 
space-time  where  appropriate  Riemann  invariants  hold,  permitting  the  integration  of  the  complete 
flow/chemistry  evolution  equations  in  an  accurate,  robust,  and  mathematicaUy  consistent  manner. 
This  method  requires  no  flux-splitting,  or  other  such  artifice,  as  is  usuaUy  required  in  conventional 
solutions  of  these  equations. 

In  order  to  allow  comparison  to  documented  results,  the  first  test  runs  with  this  unspUt  algo¬ 
rithm  were  performed  in  one  spatial  dimension  (plus  time),  for  chemicaUy-reacting  systems  where 
conventional  methods  are  known  to  perform  weU,  e.g.,  for  cases  where  the  chemical-production 


**  Case  N31N17  of  Papamoschou  (1989). 
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Fig.  IIL4.  Temporal  profile  of  dimensionless  shock  pressure  for  an  overdriven  ZND  detonation  (over¬ 
drive  factor  1.6).  Note  well-defined  single  temporal  mode. 

source  term  is  not  particularly  stiff,  and/or  cases  where  the  actual  system  is  not  characterized 
by  multi-mode  instabilities. t  The  algorithm,  however,  is  genuinely  multidimensional  and,  conse¬ 
quently,  capable  of  computing  up  to  full  four-dimensional  flows,  z.e.,  three  spatial  dimensions  and 
one  temporal  dimension.  Our  test  cases,  of  various  steady  and  unsteady  detonation  waves,  are 
in  excellent  agreement  with  well-known  results  in  the  literature,  e,g.,  Fickett  fe  Wood  1966  and 
Bourlioux  et  al  1991.  A  typical  example  of  these  test  computations  is  shown  in  Fig.  III.4.  This 
case,  an  overdriven  ZePdovich-vonNeumann-Doering  (ZND)  detonation,  with  an  overdrive  factor  of 
1.6,  is  well-known  to  be  linearly  unstable,  with  one  unstable  mode  (the  overdrive  factor  is  defined 
as  the  square  of  the  ratio  of  the  wave  speed  to  the  Chapman-Jouguet  wave  speed).  As  seen  in 
Fig.  III.4,  our  result  is  consistent  with  this  analysis.  In  this  simulation,  the  mesh  resolution  was 
set  at  20  cells  per  reaction  length. 

We  have  also  performed  simulations  of  detonations  near  the  Chapman-Jouguet  (C-J)  point,  z.e., 
with  near-unity  overdrive  factors.  These  cases  are  typically  very  difficult  to  compute,  because  of 
the  stiffness  and  other  characteristics  of  the  governing  equations  induced  by  the  terms  attributable 


We  note  that  such  systems  are  generally  linearly  and  non-linearly  unstable.  In  particular,  it  is  imperative  that 
the  instability  characteristics  of  the  numerical-simulation  model,  used  to  mimic  the  actual  system  behavior,  be 
characterized  by  as  nearly  similar  a  dynamical  system  behavior  as  the  actual  system.  These  systems,  however, 
are  also  typically  chaotic  with  a  large  number  of  (modal)  degrees  of  freedom  and  extreme  sensitivity  to  initial 
conditions  and  the  details  of  the  simulated  evolution  path.  Unfortunately,  general  proofs  for  the  requirements  to 
satisfy  this  important  criterion  are  not  available  at  this  time. 
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Fig.  IIL5.  Temporal  profile  of  dimensionless  shock  pressure  for  an  overdriven  ZND  detonation  (over¬ 
drive  factor  1.1).  Note  large  range  of  temporal  scales  present  in  the  solution. 

to  the  chemical  reactions  (c/.  footnote  and  related  discussion  above).  Consequently,  the  necessary 
resolution  often  exceeds  available  computational  resources.  The  results  of  a  stability  analysis  by  Lee 
k  Stewart  (1990)  suggest  that  in  one-dimensional  detonations  near  the  C-J  point,  we  should  expect 
multi-mode  instabilities,  i.e.,  multiple  temporal  scales,  in  contrast  to  the  single  mode  behavior 
exhibited  in  Fig.  III.4.  Our  simulations  reveal  multi-mode,  quasi-chaotic  behavior,  characterized 
by  the  creation  of  pockets  of  reaction  behind  the  leading  shock  front,  accompanied  by  the  formation 
of  secondary  shock  wave-systems  with  very  large  (local)  heat  release.  An  example  of  such  behavior 
is  presented  in  Fig.  IIL5  as  the  temporal  history  of  the  shock  pressure  for  a  ZND  detonation  with 
an  overdrive  factor  of  1.1.  For  this  case,  a  mesh  resolution  of  200  cells  per  reaction  length  was 
required  in  order  to  overcome  the  stiffness  difficulties. 

We  should  note  here  that  the  problem  of  stiffness  is  real  and  not  circumvented  by  the  time- 
explicit  algorithms  presently  available,  including  the  new  one.  This  part  of  the  problem  must  also 
be  addressed  by  the  new  algorithm  by  relying  on  adequate  time/space  resolution  (with  a  coupling 
through  the  CFL  condition).  The  improvement  here,  however,  is  realized  by  the  much  higher 
fidelity  with  which  the  system  evolution  can  be  tracked,  through  the  use  of  (the  proper)  Riemann 
invariants  in  each  case.  This  feature  of  the  new  algorithm,  at  least,  can  address  the  sensitivity  of 
the  chaotic-system  behavior  on  the  phase-space  path. 

Spatial  profiles  of  pressure  and  density,  at  a  dimensionless  time  of  /  =  10,  are  shown  in  Fig.  III. 6 
and  Fig.  III. 7,  respectively.  These  data  indicate  the  presence  of  secondary  shock  waves  and  asso¬ 
ciated  spatial  non-uniformities,  generated  by  combustion  of  previously-unreacted  material,  behind 
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Fig.  IIL6.  Spatial  profile  of  dimensionless  pressure,  at  a  dimensionless  time  of  t  =  10,  of  an  over¬ 
driven  ZND  detonation  (overdrive  factor  1.1). 

the  leading-wave  region.  Our  results  are  in  substantial  variance  and  qualitative  disagreement  with 
the  recently-published  results  by  He  &  Lee  (1995),  who  have  attempted  a  simulation  of  this  problem 
without  regard  to  the  considerations  discussed  above. 

We  are  presently  implementing  the  unsplit  algorithm  for  multidimensional  reacting  flows,  and 
expect  numerical  results  in  the  near  future. 

This  effort  is  part  of  the  Ph.D.  research  of  Miltos  Papalexandris. 

11,2.3  Image  correlation  velocimetry 

Image  Correlation  Velocimetry  (ICV)  is  a  numerical  method  which  estimates  the  velocity  field  of 
a  flow  from  a  sequence  of  images  that  are  closely  spaced  in  time.  The  method  can  analyze  particle- 
image  (PIV/DPIV)  data,  continuous-tone  images  of  convected  scalars,  or  any  image  sequence  of  a 
convected  Lagrangian  marker.  See  Tokumaru  &  Dimotakis  (1993)  for  a  more  complete  discussion 
and  references.  In  the  case  of  a  pair  of  images,  this  is  accomplished  by  iteratively  minimizing  the 
difference  of  a  mapping  of  the  first  image  and  the  second  image.  The  minimization  determines  the 
mapping  field,  ^(x),  and  is  carried  out  in  a  least-squares  sense  (equation  3a). 

min  /  {Eol^Cx)]  -  A[x]  }^d^x  ,  (3a) 


where, 


(3b) 
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Fig.  III. 7.  Spatial  profile  of  dimensionless  density,  at  a  dimensionless  time  of  t  =  10,  for  an  over¬ 
driven  ZND  detonation  (overdrive  factor  1.1). 

A  general  mapping  (displacement)  field,  ^(x),  is  described  parametrically  (equation  3b),  as  will  be 
illustrated  below. 

The  velocity  field  is  then  derived  as  the  time  derivative  of  the  displacement  field,  estimated,  for 
a  pair  of  images  at  <  =  in  turn,  as  a  first-order  time  difference,  i.e., 

U(x)  (4) 

*1  -  to 

Higher-order  estimates  in  time  are  possible  if  a  sequence,  rather  than  just  a  pair,  of  images  is 
available. 

The  details  of  the  velocity  field  representation  have  a  significant  impact  on  the  accuracy  of 
the  results  as  well  as  the  computational  efficiency.  The  current  effort  has  developed  a  spline 
representation  of  the  velocity  field,  along  with  a  new  initialization  scheme  to  increase  computational 
efficiency  and  robustness  of  the  method.  These  new  methods  have  been  used  to  successfully  recover 
the  velocity  field  of  a  simulated  Oseen  vortex  from  a  pair  of  scalar  images.  A  first  use  of  this  method 
is  also  in  progress  in  a  set  of  experiments  that  measure  the  velocity  and  vorticity  fields  generated 
by  an  impulsively-started  NACA-0012  airfoil,  at  a  high  angle  of  attack,  to  provide  feasibility  tests 
of  the  new  method.  The  latter  effort  will  be  described  in  subsequent  reports  and  publications. 

The  first  version  of  the  ICV  methodology  represented  the  displacement  (velocity)  field  by  local 
Taylor  expansions  (Tokumaru  and  Dimotakis  1993).  These  were  defined  on  a  mosaic-like  grid  over 
the  image.  This  representation  was  convenient  in  that  only  local  regions  of  the  two  images  were 
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considered  when  minimizing  the  parameters  of  the  local  displacement  expansions.  The  resulting 
velocity  field,  however,  was  not  continuous  at  the  intersections  of  adjoining  Taylor-expansion  re¬ 
gions,  An  added  term  to  minimize  discontinuities  in  the  displacement  field  and  its  lower  derivatives 
was  included  in  the  overall  cost  function  and  was  part  of  the  global  minimization. 

The  motivation  for  using  splines  in  the  current  effort  is  to  enforce  continuity  of  the  velocity 
field  rigorously  and  by  construction,  to  the  desired  order,  so  that  computational  effort  is  not 
expended  on  a  degree  of  freedom  that  will  be  removed  for  other  reasons.  Spline  choices  include 
piecewise-interpolating  polynomials.  These  are  polynomials  whose  coefficients  are  defined  over  a 
given  portion  of  the  spline  curve.  The  resulting  curve  (or  surface)  interpolates  the  data  which 
define  them.  Another  feature  of  these  splines  is  that  they  are  “global.”  Information  from  the  entire 
domain  is  necessary  to  evaluate  any  local  region.  As  with  splines  in  general,  they  are  constructed  to 
have  a  specified  degree  of  smoothness.  For  the  natural  cubic  spline,  up  to  second  order  derivatives 
are  continuous,  z.e.,  they  are  functions.  Boundary  conditions  must  be  imposed  to  evaluate 
these  splines.  In  the  case  of  the  natural  cubic  spline,  the  procedure  is  closed  by  setting  the  second 
derivatives  to  zero  at  the  boundaries. 

While  piecewise-interpolating  polynomials  have  proven  useful  in  many  applications,  they  have 
serious  shortcomings  for  use  in  ICV.  The  minimization  scheme  of  ICV  relies  on  repeated  evalu¬ 
ations  of  the  difference  between  the  mapped  first  image  and  the  second  image,  as  the  mapping 
(displacement /velocity)  field  parameters  are  varied.  With  global  splines,  information  from  the  en¬ 
tire  domain  is  required  to  evaluate  even  a  small  region.  The  many  matrix  inversions  necessary  to 
solve  for  the  coefficients  of  the  spline  result  in  a  substantial  computational  burden.  Furthermore, 
splines  which  interpolate,  z.e.,  go  through  the  actual  data  (control)  points,  are  generally  prone  to 
large  undulations  between  the  data.  This  can  in  turn,  lead  to  instabilities  in  the  ICV  algorithm. 

These  two  problems  are  addressed  in  the  current  ICV  implementation  through  the  use  of  B- 
splines  (c/.  Bartels  et  aL  1987),  These  splines  only  approximate  the  data  which  define  them. 
The  relaxation  of  the  requirement  of  strict  interpolation  eliminates  the  “overshoot/undershoot” 
tendencies  of  the  piecewise  interpolating  natural  bicubic  splines  discussed  above.  Furthermore, 
B-splines  are  “local,”  «.e.,  only  a  small  neighborhood  of  surrounding  data  is  necessary  to  evaluate  a 
local  region.  This  feature  results  in  substantial  computational  savings,  especially  when  considering 
surfaces  defined  in  terms  of  many  points. 

An  important  question  that  arises  is  what  order  spline  is  required  to  represent  an  actual  velocity 
field.  Consider  the  Navier-Stokes  momentum  and  vorticity  equations,  z.e.. 


du  du 

duj  du: 


1  dp  _2 
-  -  ^  +  Z/V^u 
p  ox 


(-4 


u  +  z/V^u?  , 


(5) 


for  incompressible  flow.  As  can  be  seen  from  the  vorticity  equation,  which  imposes  the  stricter 
continuity  requirements  of  the  two,  the  vorticity  field,  u?(x),  must  be  at  least  C^(x)-continuous.  In 
fact,  a  C^(x)  u?-field  would  quickly  acquire  continuous,  higher-order,  spatial  derivatives.  To  limit 
the  continuity  constraints  to  a  manageable  order  we  have  elected  to  represent  (incompressible) 
u-fields  as  C^(x),  resulting  in  vorticity  fields  that  are  C^(x).  This  provides  a  two-order  continuity 
margin  with  respect  to  the  momentum  equation  and  one-order  continuity  margin  with  respect  to 
the  vorticity  equation. 
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What  follows  is  a  brief  description  of  univariate  B-splines,  i.e.,  for  curves  that  are  functions  of 
one  independent  variable.  The  univariate  case  is  easily  generalized  into  the  bivariate  case,  as  used 
in  ICV,  i.e.,  surfaces  that  are  functions  of  two  independent  variables,  say,  of  x  and  y,  as  appropriate 
for  two-dimensional  fields.  For  a  more  complete  description  of  B-splines,  see  (Bartels  et  al.  1987) 


Fig.  III.8.  Equally-weighted  basis  functions  of  the  quintic  B-spUne. 


Fig.  III. 9.  Simple  curve  generated  with  quintic  B-splines. 


In  the  most  general  sense,  a  B-spline  can  be  thought  of  as  a  summation  of  basis  functions  over 
a  given  region.  Figures  III. 8  and  equationIII.9  show  the  basis  functions  for  the  quintic-order  spline 
as  the  beU-curve-shaped  solid  lines  at  the  bottom.  Basis  functions  are  identically  zero  everywhere 
(compact  support),  except  for  the  local  region  centered  around  their  corresponding  data  (control) 
point.  As  illustrated  in  Fig.  III. 9,  the  basis  functions  are  weighted  by  their  corresponding  data 
then  summed  over  a  region.  The  result  of  the  summation  is  the  desired  spline  curve,  shown  as  the 
broken  lines  approximating  the  control  points  above  the  basis  functions. 

One  should  also  note  that  the  smoothness  of  the  resulting  curve  is  independent  of  the  irregularity 
of  the  data.  Furthermore,  the  curve  will  always  lie  in  a  region  bounded  by  the  data  points.  The 
ability  of  B-splines  to  produce  smooth  curves  is  a  direct  consequence  of  how  the  basis  functions 
are  constructed.  While  there  are  many  ways  to  derive  the  B-spline  basis  functions,  the  following 
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recursion  formula  is  a  concise  mathematical  definition,  for  functions  centered  at  the  origin. 

^0/^^  ^  5  for  -1/2  <  a;  <  1/2  , 

^  ^  1 0 ,  otherwise; 

5”(a;)  =  , 

where  is  the  n*^^-order  basis  function  and  denotes  convolution.  The  l®*'-order  function 

is  found  recursively  by  convolving  the  0''*^-order  function  with  itself.  The  2"^-order  is  found  by 
convolving  the  0*'**-order  with  the  l®*’-order,  etc.  This  produces  basis  functions  that  are  centered  at 
a;  zz  0,  which  may  be  translated  in  x  to  position  them  as  necessary.  The  spline  curve,  S{x)^  is  then 
given  by, 

k 

S(x)  =  Y,q,Bf(x)  ,  (7) 

z=0 

where  the  qi  are  the  weights  of  the  control  points  and  define  the  curve,  parametrically,  for  a 
given  set  of  basis  functions. 

In  Figs.  IIL8  and  III. 9,  the  curve  is  only  shown  for  the  inner  regions  of  the  control  points,  i.e., 
between  5  and  7  on  the  independent- variable  axis.  Quintic  splines  require  the  support  of  six  basis 
functions  to  be  completely  defined.  Various  schemes  exist  to  accommodate  the  limited  support  at 
the  boundaries,  which  wiU  not  be  covered  here.  One  such  scheme  has  been  used  in  the  current  ICV 
method,  allowing  the  surface  to  extend  to  the  boundary  data  and/or  accommodate  appropriate 
boundary  conditions,  in  each  case. 


Fig.  III.  10.  Single  basis  function  of  biquintic  B-sphne. 

Generalization  to  two-dimensions  can  be  realized  using  biquintic  B-splines.  The  biquintic  basis 
function  in  Fig.  III. 10  was  defined  by  the  tensor  product  of  one-dimensional  quintic  basis  (equa¬ 
tion  6).  Some  lower-order  tensor-product-defined  bivariate  B-spline  bases  are  given  in  equation  8. 
The  basis  is  the  two-dimensional  top-hat  function  and  interpolation  using  B-splines  of  this 
order  results  in  a  ^^nearest  neighbor”  approximation. 


1,  for-l/2<x,y<l/2; 

0,  otherwise. 


(8a) 
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Bilinear  B-splines  interpolate  four  data  points  in  a  bilinear-sense.  The  equation  for  the  basis 
function  is  given  explicitly  as  foDows, 

'  (1  +  a:)(l  +  y)  for  -1  <  x,y 

(1  -  a;)(l  +  y)  for  0  <  <  1,  -1  <  ^  <  0, 

jB  '  (x,  y)  =  {  (1  —  5;)(1  ~  y)  for  0  <  rr,  2/  <  1,  (8b) 

(1  +  a:)(l  —  y)  for  — 1  <  a;  <  0,  0  <  2/  <  1, 

0  otherwise. 


Fig.  III. 11.  Laser-induced  fluorescence,  512  x  512-image  data  used  for  Oseen  vortex  test  case.  Data 
digitized  at  llbits/pixel. 

The  new  ICV  algorithm  has  been  tested  on  a  the  model  flow  field  of  an  Oseen  vortex,  which 
is  an  analytical  solution  for  the  temporal  decay  of  a  vortex  filament.  The  analytical  velocity  field 
was  used  to  map  the  image  in  Fig.  III.  11  into  a  simulated  scalar  image  that  would  occur  at  a  later 
time  if  the  first-image  scalar  field  was  convected  by  the  corresponding,  two-dimensional  field.  The 
512  X  512  test  image  in  Fig.  III.ll  was  formed  using  a  fluorescent  dye  and  a  laser  sheet  in  water, 
digitized  at  llbits/pixel.  These  two  images  were  then  used  as  the  input  to  the  ICV  algorithm. 
The  Oseen  vortex,  defined  in  equation  9,  was  centered  in  the  image, 

Urir^t)  =  0  , 


(9) 
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with  the  parameters  C  and  ft  chosen  to  produce  an  interesting  velocity  distribution.  Fignres  III.  12 
and  III. 13  are  surface  plots  of  the  vertical  component  of  velocity,  v(x,y),  and  vorticity,  oj{x,y), 
respectively. 


Fig.  III. 12.  Surface  plot  of  v{x,y)  for  Oseen  vortex  used  to  map  test  image. 

The  ICV  algorithm  iteratively  minimizes  the  difference  between  the  mapped  first  image  and 
the  second  image  (equation  3a).  A  new  initialization  routine,  based  on  Fourier  cross-correlation 
methods,  similar  to  DPIV  techniques  {e.g.,  Willert  and  Gharib  1991),  was  developed  to  generate  a 
first  estimate  of  the  velocity  field,  which  is  then  used  as  the  input  to  the  main  minimization  routine. 
For  DPIV  (particle)  images,  cross-correlation  region  window  sizes  are  of  the  order  of  32  X  32  pixels. 
For  the  purpose  of  correlating  images  of  scalars,  window  sizes  on  the  order  of  64  X  64,  to  128  X  128, 
were  required  for  a  reliable  initial  estimation. 

The  new  initialization  routine  goes  a  step  beyond  standard  DPIV  methods  in  that  it  generates 
a  smooth  displacement  field  by  construction  and  is  an  iterative  scheme.  The  velocity  field  from  the 
initial  Fourier  cross  correlation  is  represented  using  biquintic  B-splines  and  then  used  to  produce  a 
mapping  of  the  first  image.  One  can  iterate  for  a  more  accurate  solution  by  performing  the  cross 
correlation  again.  In  this  iteration,  the  mapping  of  the  first  image  is  correlated  with  the  second 
image.  This  second  cross-correlation  iteration  is  able  to  solve  for  translational  errors  in  the  initial 
mapping.  In  this  manner,  one  arrives  at  a  correction  that  is  added  to  the  previous  mapping. 

The  initial  estimation  of  the  velocity  field  is  then  passed  to  the  main  ICV  routine.  The  scheme  is 
described  in  detail  by  Tokumarn  and  Dimotakis  (1993),  except  that,  in  the  present  implementation, 
B-splines  are  used  instead  of  local  Taylor-expansions.  Modifications  in  the  implementation  of  the 
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Fig.  111.13.  Surface  plot  of  vorticity,  uj{x,y)  for  Oseen  vortex  used  to  map  test  image. 


Fig.  III.14.  Surface  plot  of  the  vertical  component  of  velocity,  v{x,y),  deduced  from  ICV  solution 
for  the  simulated  Oseen  vortex. 

original  algorithm  have  been  performed  to  take  advantage  of  the  B-spline  velocity-field  representa¬ 
tion.  In  particular,  for  a  given  region  of  the  two  images  to  be  correlated,  the  minimization  routine 
is  not  completely  local  in  that  it  is  allowed  to  vary  the  neighboring  control  points  of  the  B-spline 
representation  of  the  velocity  field  to  achieve  the  best  mapping.  Results  of  the  ICV  minimization 
are  shown  in  Fig.  III.14  and  Fig.  III.15. 
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Fig.  IIL15.  Surface  plot  of  vorticity,  uj{x^y)^  deduced  from  ICV  solution  for  the  simulated  Oseen 
vortex. 

The  error  surface  for  the  experimentally-determined  Uexp ( 2:,  y)- velocity  field,  z.e., 


^v{x,y) 


-  Vthix,y) 

max{|?;(a:,j/)|} 


(10) 


is  plotted  in  Fig.  III. 16,  expressed  as  a  percent  of  the  maximum  u-velocity  in  the  field,  f.e.,  = 

100  6^,(3;,  ?/).  The  error  in  the  deduced  u- velocity  is  0.4%  (rms)  over  the  whole  field  and  everywhere 
less  than  1.2%. 


Fig.  III.17.  Surface  plot  of  percent  error  in  ICV  solution  for  the  vorticity,  u{x,y),  of  the  simulated 
Oseen  vortex  (c/.  equation  11  and  related  discussion). 
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Fig.  IIL16. 


Surface  plot  of  percent  error  in  ICV  solution  for  the  vertical  velocity,  v{x^y)^  of  the 
simulated  Oseen  vortex  (c/.  equation  10  and  related  discussion). 


The  error  surface  for  the  experimentally- determined  vorticity  field,  2/)? 


eui{x,y) 


^exp(a:,y) 

max{|u;(x,?/)|} 


(11) 


is  similarly  plotted  in  Fig.  111.17,  as  a  percent  of  the  peak  vorticity  in  the  flow  field,  2.e.,  2r  = 
100£a;(^5  2/)-  The  error  in  the  deduced  vorticity  over  the  whole  field  is  0.8%  (rms)  and  everywhere 
less  than  2.1%,  for  this  test  case,  with  much  of  it,  as  can  be  seen,  traceable  to  the  boundaries  of 
the  ICV  flow  domain. 

In  work  in  progress,  we  are  processing  the  data  from  the  experiment  investigating  the  flow  over 
the  impulsively-started  airfoil.  This  part  of  the  eifort  wiU  require  further  development  of  the  B- 
spline  flow-field  surface  description  and  the  ICV  method  to  allow  more  complex  geometries  as  weU 
as  flow/body  boundary  conditions  to  be  accommodated. 

The  important  extension  of  this  work  to  compressible  flow  fields  that  are  not  necessarily  smooth, 
e.g.^  flows  with  shocks,  is  also  part  of  the  overall  efibrt.  Envisaged  algorithms  will  attempt  to  fit 
shocks  first  and  then  solve  for  piecewise-smooth  velocity-field  regions,  as  described  above. 

This  effort  is  part  of  the  Ph.D.  research  of  Galen  Gornowicz. 
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CHAPTER  IV 

CHEMISTRY  IN  NONUNIFORM  FLOW 

IV.l  DISSOCIATION  RATES  WITH  VIBRATIONAL  NONEQUILIBRIUM 

Objectives  and  Status  of  Research 

One  of  the  primary  uncertainties  in  the  modeling  of  hypersonic  reacting  flows  is  how  the  vibra¬ 
tional  state  of  a  molecule  affects  its  dissociation  rate.  Therefore,  the  main  objective  of  this  research 
is  to  test  chemical  reaction  rate  models  in  the  presence  vibrational  nonequilibrium.  Our  initial  work 
showed  that  the  flow  of  nitrogen  over  spheres  at  typical  T5  conditions  is  not  very  sensitive  to  the 
choice  of  several  popular  vibration-dissociation  coupling  models.  During  the  last  year,  we  have  been 
using  computational  fluid  dynamics  to  design  a  new  experiment  that  shows  strong  sensitivity  to  the 
coupling  model.  We  have  found  several  geometries  and  flow  conditions  that  are  good  candidates; 
they  wiU  be  tested  in  T5  during  the  summer.  This  experiment  will  provide  the  first  data  that  can 
be  used  for  the  direct  validation  of  vibration-dissociation  coupling  models,  and  it  wiU  provide  new 
insight  into  how  the  vibrational  state  affects  the  dissociation  rate. 

The  new  experimental  configuration  is  a  double-wedge  or  double-cone  geometry;  the  first  cone 
or  wedge  causes  an  attached  shock  wave  to  form,  and  interact  with  the  detached  shock  from  the 
second  wedge  or  cone.  This  creates  a  complicated  flow  field  that  is  very  sensitive  to  the  vibration- 
dissociation  coupling  model.  For  example.  Fig.  IV.l  plots  the  pressure  contours  in  the  flow  field 
of  a  double- cone  geometry  at  typical  T5  conditions.  The  shock  wave-boundary  layer  interaction 
produces  a  region  of  separated  flow  at  the  intersection  of  the  cones.  Also,  the  shock-shock  interaction 
produces  a  dramatic  pressure  rise,  followed  by  an  oscillating  pressure  distribution  on  the  surface  of 
the  second  cone.  The  only  difference  between  these  two  calculations  is  the  vibration-dissociation 
coupling  model.  The  Park  coupling  model  produces  a  flow  field  with  a  much  smaller  separation 
zone  than  the  CVDV  model.  The  difference  in  the  separation  zone  causes  a  large  change  in  the 
surface  pressure  distribution,  as  seen  in  Fig.  IV.2.  There  is  a  similar  change  in  the  convective 
heating  rate  distribution.  The  experimental  models  will  be  instrumented  with  pressure  and  heat 
transfer  gauges,  and  Mach-Zehnder  interferograms  will  be  made. 

During  the  design  of  the  double-wedge/cone  experiments,  we  discovered  some  interesting  flow 
features  and  what  appears  to  be  a  new  shock-shock  interaction.  For  example,  the  Mach  6.9  flow  of 
a  perfect  gas  over  a  15°-58°  double-wedge  creates  a  shock-shock  interaction  that  is  similar  to,  but 
different  than,  a  Type  IV  interaction  (Compare  with  the  shock-shock  interaction  of  section  II.2).  As 
seen  in  Fig.  IV. 3,  the  shock  that  is  transmitted  from  the  triple  point  undergoes  a  regular  reflection 
from  the  first  wedge  surface,  creating  a  supersonic  flow  along  the  surface  of  the  body.  This  flow  is 
similar  to  an  underexpanded  jet,  making  it  undergo  alternating  expansions  and  recompressions. 


(a)  Park  \/T%  coupling  model.  (b)  CVDV  coupling  model. 

Figure  IV.  1.  Pressure  contours  for  double-cone  geometry  at  typical  re-entry  conditions  (p 
0.01  kg/m^,  Too  =  2260  K,  Uoo  =  6690  m/s,  base  height  is  2 inches  (5.08  cm)). 


Distance  along  Body  (mm) 

Figure  IV.2.  Surface  pressure  distribution  on  double-cone  geometry  shown  in  Fig.  IV. 1. 
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Figure  IV. 3.  Mach  number  and  pressure  contours  of  the  new  shock  interaction.  The  flow  is  inviscid 
with  Moo  =  6.9  and  7  =  1.4;  computation  performed  on  a  512  x  512  grid. 
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IV.2  ELECTRON-DRIVEN  REACTIONS  IN  HYPERSONIC  FLOW 

Background  and  Objective 

Electrons  resulting  from  ionization  due  to  bow-shock  heating  play  an  important  role  in  deter¬ 
mining  the  physical  and  chemical  properties  of  hypersonic  flows. 

Processes  of  interest  include  vibrational  excitation,  formation  of  negative  ions  (which  may  then 
fragment),  and  electron  impact  excitation  and  deexcitation  of  excited  metastable  states.  Further¬ 
more,  with  temperatures  as  high  as  10,000  Kelvin  due  to  non-equilibrium  overshoots,  electrons  in 
the  tail  of  the  energy  distribution  function  can  lead  to  electronic  excitation  of  ground  state  species, 
particularly  those  that  are  vibrationally  excited.  The  available  data  base  for  these  electron-molecule 
collision  processes  which  can  occur  in  flowfields  is  generally  fragmentary,  and,  in  many  cases,  in¬ 
consistent.  Furthermore,  cross  sections  for  electronic  excitation  of  molecules  by  electrons  with 
barely  enough  energy  to  drive  the  excitation  and  for  electron  collisions  with  electronically  excited 
molecules  and  with  radicals  such  as  OH  are  essentially  non-existent  and  experimentally  inaccessible. 
Finally,  the  experimental  effort  to  measure  electron-molecule  collision  cross  sections  is  very  limited 
and  shrinking  in  the  United  States. 

An  obvious  strategy  for  obtaining  these  electron  collision  cross  sections  is  to  calculate  them. 
However,  the  significant  computational  demands  of  such  calculations  have  made  such  an  approach 
impractical  to  date  with  conventional  (sequential)  computers.  The  objective  of  our  work  is  to 
exploit  the  high-performance  computing  provided  by  parallel  computers,  consisting  of  arrays  of  com¬ 
mercially  available  microprocessors,  to  calculate  the  cross  sections  for  electron-molecule  collisions 
needed  in  modelling  the  chemical  and  physical  properties  of  hypervelocity  flows.  With  appropriate 
computational  strategies,  the  high  aggregate  speeds  and  large  memory  of  these  parallel  computers 
make  such  an  approach  feasible. 


Progress  and  Highlights 

At  the  relatively  low  impact  energies  of  interest  in  hypersonic  flowfields,  the  accurate  treatment 
of  the  electron-molecule  collision  process  is  a  computationally  challenging  problem.  In  these  studies 
we  employ  a  multichannel  extension  of  the  variational  principle  for  scattering  amplitudes  originally 
introduced  by  Schwinger.^  This  formulation,  known  as  the  Schwinger  multichannel  method,^  was 
specifically  designed  to  address  low-energy  electron-molecule  collisions. 

As  in  most  variational  approaches,  there  are  two  major  computational  tasks  in  these  studies:  (i) 
calculation  of  the  primitive  matrix  elements  of  the  operators  arising  in  the  variational  expression 
for  a  given  expansion  basis,  and  (n)  transformation  of  these  primitive  matrix  elements  to  gener¬ 
ate  the  physical  matrix  elements  of  the  linear  equations  defining  the  variational  solution.  In  the 
Schwinger  multichannel  method,  as  in  the  original  Schwinger  method,  the  trial  wave  function  need 
not  satisfy  scattering  boundary  conditions;  square-integrable  functions  such  as  Cartesian  Gaussians 
may  be  employed.  In  this  basis  all  matrix  elements  in  the  variational  expression  can  be  evaluated 
analytically,  except  those  involving  the  Green’s  function.  These  Green’s  function  matrix  elements 
have  no  known  analytic  form  and  must  be  done  by  quadrature.  In  fact,  evaluation  of  these  matrix 
elements  becomes  the  computationally  intensive  step  of  our  calculations.  The  occurrence  of  these 
difficult  matrix  elements  in  variational  principles  based  on  integral  equations  has  discouraged  their 
use  to  date,  in  spite  of  highly  desirable  features  of  these  approaches  such  as  the  flexibility  of  em- 
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ploying  trial  functions  which  do  not  satisfy  scattering  boundary  conditions  and  their  very  favorable 
convergence. 

Though  initially  implemented  on  conventional  (sequential)  computers,  the  computationally  in¬ 
tensive  step  of  these  calculations  is  ideally  suited  to  computers  of  distributed-memory  parallel 
architectures.  Its  most  compute-intensive  tasks  are  (i)  the  evaluation  of  a  very  large  number  (typ¬ 
ically  or  more)  of  elementary  integrals  describing  the  interaction  of  pairs  of  electrons,  a  step 
requiring  a  code  of  about  2000  lines  of  Fortran,  and  {ii)  the  transformation  of  these  elementary 
quantities  into  the  matrix  elements  appearing  in  the  variational  expression.  The  evaluation  of  the 
two-electron  integrals  is  readily  parallelized:  one  simply  has  each  processor  independently  calculate 
a  different  subset  of  the  integrals.  Step  (n),  the  transformation  of  the  integrals,  can  be  organized 
as  a  series  of  large  distributed  matrix  multiplications  between  the  global  array  of  integrals  and  a 

o 

distributed  transformation  matrix  that  is  easy  to  construct.'^ 

The  overall  performance  of  our  program  on  the  Intel  Paragon  (512  i860  processor)  or  CRAY 
T3D  (256  DEC  alphas)  cannot  be  characterized  by  a  single  number.  Although  evaluating  and 
subsequently  combining  the  integrals  are  by  far  the  dominant  steps  in  the  calculation,  the  relative 
importance  of  these  two  steps  is  strongly  dependent  on  the  specific  calculation.  Depending  on  the 
nature  of  the  calculation  at  hand,  then,  and  taking  into  account  load  balance,  communication,  I/O, 
and  other  overhead,  we  typically  see  throughout  speeds  in  the  range  from  3  to  10  GFLOP.  To  put 
these  numbers  in  perspective,  our  original  sequential  program,  which  is  only  partially  vectorizable, 
averages  about  30  MFLOP  on  a  CRAY  Y-MP  processor.  This  increase  of  over  two  orders-of- 
magnitude  in  computational  performance  significantly  enhances  our  ability  to  fill  the  large  gaps  in 
the  electron  collision  data  base  needed  for  modelling  reactive  flows  of  interest. 

We  have  exploited  this  computational  procedure  and  parallel  computer  resources  at  the  Institute 
and  Wright-Patterson  and  Eglin  Air  Force  Bases  to  carry  out  the  following  calculations: 

•  We  have  studied  the  cross  sections  for  near-threshold  excitation  of  the  lowest  lying  excited 
state  of  OH  (A^S^)  by  low-energy  electrons.  OH  fragments  are  produced  in  flows  of  interest  via 
dissociation  of  water  vapor.  No  measured  data  is  available  for  electron  impact  on  the  OH  fragment. 

•  We  have  completed  calculations  of  the  cross  sections  for  electron  impact  excitation  of  the  A 
state  of  NO.  This  upper  state  is  responsible  for  the  7-band  emission  which  is  prominent  in  the 

bow  shock  of  flows  over  hardbodies.  There  have  also  been  no  measurements  of  these  cross  sections. 

•  We  have  completed  extensive  calculations  of  the  cross  sections  for  electron  impact  excitation 
of  five  excited  electronic  states  of  CO2.  There  is  again  essentially  no  measured  data  for  any  of  these 
cross  sections. 


References 

1.  J.  Schwinger,  Phys.  Rev.  72,  742  (1947) 

2.  K.  Takatsuka  and  V.  McKoy,  Phys.  Rev.  A  30,  1734  (1984);  C.  Winstead  and  V.  McKoy,  Phys. 
Rev.  A  5,  648  (1972) 

3.  C.  Winstead,  P.G.  Hipes,  M.  A.  P.  Lima,  and  V.  McKoy,  J.  Chem.  Phys.  94,  5455  (1991) 


-  60- 


Publications 

The  following  publications  have  resulted  from  the  research  outlined  above: 

1.  Lee,  C-H,  Winstead,  C.,  Pritchard,  H.P.,  and  McKoy,  V.  Studies  of  Collisions  of  Low-Energy 
Electrons  with  CO2?  Phys  Rev.  A  (submitted). 

2.  Zobel,  J.,  Mayer,  U.,  Jung,  K.,  Ehrhardt,  H.,  and  Pritchard,  H.,  Winstead,  C.,  and  McKoy, 
V.,  Absolute  Cross  Sections  for  Electron- Impact  Excitation  of  CO  Near-Threshold^  J.  Phys.  B. 
(submitted). 

3.  Winstead,  C.  and  McKoy,  V.,  Studies  of  Electron- Molecule  Collisions  on  Massively  Parallel 
Computers^  Modern  Electronic  Structure  Theory,  Edited  by  D.  Yarkony  (World  Scientific,  1995). 

4.  Winstead,  C.,  Pritchard,  H.P.,  and  McKoy,  V.,  Electron- Molecule  Collisions  for  Plasma  Mod¬ 
elling  Experiences  on  the  CRAY  T3D  CRAY  CHANNELS  (Cray  Research,  Inc.,  1995). 

5.  Winstead,  C.,  Pritchard,  H.P.,  and  McKoy,  V.,  Parallel  Computations  of  Electron- Molecule 
Collisions,  Invited  Article  in  Computational  Science  and  Engineering,  IEEE  Computer  Society 
(1995). 


Personnel  Associated  with  Research 

Principal  Investigator:  Vincent  McKoy,  Professor  of  Theoretical  Chemistry 

One  graduate  student  (C-H.  Lee)  and  a  Senior  Research  Fellow  (C.  Winstead,  part-time)  are 
involved  in  this  effort. 


Presentations 

The  research  outlined  above  has  been  presented  in  the  following  talks: 

•  Electron- Molecule  Collisions  for  Modelling  Plasmas,  Invited  Talk  at  the  Mardi  Gras  ’94  Con¬ 
ference  “Toward  Teraflop  Computing  and  New  Grand  Challenge  Applications”,  Baton  Rouge,  1994. 

•  Electron- Molecule  Collisions  for  Modelling  Plasmas,  Invited  Talk  at  the  High-Performance 
Computing  Symposium,  1994,  Simulation  MultiConference,  La  JoUa,  1994. 

•  Studies  of  Electron- Molecule  Collisions  on  Highly  Parallel  Computers,  Joint  Chemical  Physics- 
Mechanical  Engineering  Seminar,  Ohio  State  University,  May  1994. 

•  Data  for  Modelling  Materials-Processing  Plasmas:  The  Impact  of  Parallel  Computers,  Invited 
Talk  at  the  Defense  Science  Study  Group  Symposium  on  Applications  of  Advanced  and  Innova¬ 
tive  Computational  Methods  to  Defense  Science  and  Engineering,  Institute  for  Defense  Analyses, 
October  1994. 

•  Data  for  Modelling  Materials-Processing  Plasmas:  The  Impact  of  Parallel  Computers,  Invited 
Talk  in  the  High-Performance  Computing  Seminar  Series,  National  Institute  for  Standards  and 
Technology,  April  1995. 


61- 


•  Studies  of  Electron- Molecule  Collisions  on  Highly  Parallel  Computers^  Bucy  Lecture,  Texas 
Tech  University,  February  1995. 

•  Studies  of  Electron- Molecule  Collisions  on  Parallel  Computers^  Invited  Talk  at  the  1995  Annual 
Meeting  of  the  Division  of  Atomic,  Molecular  and  Optical  Physics,  American  Physical  Society,  May 
1995,  Toronto,  Canada. 


-62- 


IV.3  NONEQUILIBRIUM  LEEWARD  SHOCK- VORTEX  AERODYNAMICS 


Objectives  and  Status  of  Research 

The  objectives  of  this  research  are  to  study  the  hypervelocity  flow  at  enthalpy  sufficient  to 
activate  dissociation/recombination  chemistry  about  cones  and  blunt-nosed  delta- wings  at  large 
incidence,  with  particular  reference  to  the  effect  of  chemistry  on  the  structure  of  the  leeward  flow 
shock/ vortex  system.  Significant  progress  has  been  made  in  two  relevant  areas.  These  are 

(1)  A  new  kinetic  scheme  based  on  the  equilibrium  flux  method  {EFM)  and  modified  using  Osher 
intermediate  states  has  been  developed.  This  new  scheme,  called  EFMO^  is  designed  to  retain 
the  robustness  of  the  original  Boltzmann-type  method  in  capturing  intense  shock  waves  and  at 
the  same  time  substantially  reduce  the  numerical  diffusivity  of  the  original  EFM  method  for 
resolution  of  boundary  layers. 

(2)  A  fully  implicit  version  of  the  viscous  three-dimensional  (3D)  flow  code,  PGP3D  has  been  made 
developed  with  EFMO  incorporated.  The  code  has  been  tested  extensively  on  viscous  laminar 
cone  flow  at  zero  and  at  finite  incidence. 

(3)  The  Light  hill- Freeman  ideal-dissociating  gas  chemistry  model  has  been  incorporated  into  both 
the  explicit  and  implicit  versions  of  PGP3D  and  is  currently  undergoing  testing. 

A  comprehensive  statement  of  accomplishments  in  the  areas  (1-3)  above  is  outlined  below. 
These  statements  constitute  the  current  status  of  the  research  effort.  They  are  followed  by  a  brief 
statement  of  the  expected  progress  within  the  next  year. 

IV.3.1  The  EFMO  Scheme 

The  EFMO  scheme  (EFM  and  Osher)  is  based  on  a  formal  similarity  between  the  mathematical 
expressions  of  the  EFM  numerical  flux  and  the  Osher  numerical  flux.  In  Osher’s  method  (Osher  & 
Solomon  1982),  the  numerical  flux  is  given  by 

Fo{Ul.Ur)^^I2{fl^Fr-  [  \A\dw-  [  \A\dw),  (1) 

V  Jln  Jnl  ) 

where  LN  refers  to  the  linearly  degenerated  subpath  which  corresponds  physically  to  a  contact 
discontinuity,  and  N L  refers  to  the  genuinely  nonlinear  subpaths  which  correspond  to  either  com¬ 
pression  or  rarefaction  waves.  The  matrix  |A|  is  the  absolute  value  of  the  Jacobian  matrix  of  the 
exact  flux  F  and  is  defined  in  the  following  sense: 

Ml  =  R\A\R-^  |A|  =  Diag(|Ap|)p=i,...^,  (2) 

where  Xp  is  an  eigenvalue  of  A. 

The  T'TM  numerical  flux  (PuUin  1980)  defined  traditionally  as 

Fefm{Ul,  Ur)  =  F+{Ul)  +  F-(Ur),  (3) 

can  be  put  under  a  form  similar  to  Osher’s  numerical  flux 

Fefm{Ui,Ur)=  '^/2(fi  + Fr-  f  \B\dw-  I  \B\dw\ , 

V  JLN  JNL  / 


(4) 
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where  \B\  is  now  defined  as  \B\  =  B+  -  B  where  is  the  Jacobian  matrix  of  the  fluxes  F^. 

The  idea  of  EFMO  is  to  observe  that  for  the  nonlinear  paths  (i.e.  shock  waves,  in  particular) 
EFM  performs  extremely  well  while  it  tends  to  artificially  broad  contact  discontinuities  owing 
to  its  large  numerical  diffusivity.  On  the  other  hand,  Osher  scheme  is  known  to  present  some 
inconsistencies  in  the  vicinity  of  sonic  lines  or  of  strong  shocks  while  giving  excellent  resolution  of 
contact  discontinuities.  The  present  idea  is  to  retains  that  part  of  the  integral  corresponding  to  the 
nonlinear  path  (shock  wave)  in  the  EFM  numerical  flux  while  replacing  the  integral  corresponding 
to  the  linear  path  (contact  discontinuity)  by  the  expression  from  Osher’s  numerical  flux.  The 
resulting  numerical  flux  is  then 

FeFMo{Ul,  Ur)  =  1/2  (j'l  +Fr-  \A\dw  -  \B\dw^  ,  (5) 

or  equivalently 

UeFMo{Ul-,Ur)=  Fefm{Ul^Ur)F^I2  (  {\B\-\A\)dw,  (6) 

Jln 

where  the  integral  represents  an  anti- diffusive  flux  which  almost  vanishes  in  the  vicinity  of  shocks 
and  removes  the  excessive  numerical  dissipation  near  contact  discontinuities.  In  practice,  the 
previous  relation  can  be  further  simplified  by  defining  two  intermediate  states  Ui  and  U2  defined 
by  using  Riemann  invariants  in  the  sense  of  Lax  as  implied  in  the  original  Osher’s  method. 

These  two  intermediate  states  are  then  used  to  calculate  the  anti-diffusive  flux  so  that  eventually 
the  EFMO  numerical  flux  reads 

FefMo{Vl,Ur)  =  FErM(UL,VR)  +  {  +  f (J))* '  J  J  ,  (7) 

where  is  the  common  value  of  the  normal  velocity  component  for  both  intermediate  states  and 
are  the  regular  T’FM  fluxes  which  appear  in  Eq.  (3). 

IV.3.2  PGP3D  and  Implicit  EFMO 

An  implicit  version  of  PGP3D  incorporating  EFMO  {PGP3D.2)  has  been  developed  based  on 
a  linearization  of  the  fluxes  using  an  approximate  analytical  evaluation  of  the  Jacobian  in  which 
only  the  regular  EFM  part  of  the  numerical  flux  is  retained.  This  technique  provided  a  rate 
of  convergence  with  EFMO  comparable  to  that  obtained  with  EFM,  for  which  the  Jacobian  is 
exact.  For  three-dimensional  flow  this  procedure  results  in  a  block  hepta-diagonal  system  of  linear 
equations  which  must  be  solved  at  each  time-step,  and  this  is  done  by  factorization  into  tridiagonal 
systems  which  are  solved  by  standard  methods.  The  extension  of  EFMO  to  high  order  of  accuracy  is 
achieved  using  the  usual  MUSCL  technique  which  consists  of  reconstructing  the  primitive  variables 
on  each  side  of  a  given  interface  where  the  flux  is  calculated. 

The  PGSD/FFMO  code  has  been  subject  to  extensive  testing  for  a  variety  of  supersonic  and 
hypervelocity  frozen  flows.  These  include  viscous  flow  over  a  flat  plate  at  zero  incidence  at  at  various 
Mach  numbers,  the  flow  over  a  flat-plate  ramp  combination  at  M  =  14.2  and  at  a  ramp  angle  of  24°, 
and  hypersonic  viscous  cone-flow  at  zero  and  at  finite  incidence.  The  aim  is  to  produce  a  reliable 
method  which  will  accurately  reproduce  compressible  laminar  boundary-layer  behavior  and  strong 


Figure  IV.4.  Flow  property  profiles  ?iX  xj L  =  0.1  on  a  cone  of  half-angle  10^  at  zero  incidence  in  a 
hypersonic  stream  at  M  =  7.95.  Frozen  flow,  =  4.2  X  10^  based  on  distance  L  from 
apex  and  Pr  =  1.  Adiabatic  wall,  (a)  u/Vqo  versus  rj  =  y/L  where  y  is  normal  to  the  wall, 
(b)  pIpoo  versus  p. 


(a)  (b) 


Figure  IV. 5.  For  key,  see  Figure  IV.4.  (a)  T/Too  versus  p  =  y/L,  (b)  Detail  of  (a).  Diamond  gives  wall 
temperature,  using  recovery  factor  r  =  1. 

detached  or  attached  shocks.  Figures  IV.4  and  IV.5  show  profiles  of  tangential  velocity  ufVoo^ 
static  pressure  pjpoo  and  temperature  profiles  T/Too,  normahzed  against  free-stream  conditions,  at 
a  station  xjL  =  0.1  along  a  cone  of  half- angle  10^  at  zero  incidence  in  hypersonic  viscous  laminar 
frozen  flow  at  Moo  =  7.95,  where  x  is  measured  along  the  cone  axis  and  L  is  the  cone  length. 
The  Reynolds  number  based  on  free-stream  conditions  and  L  was  4.2  x  10^.  The  Prandtl  number 
was  Pr  =  1  and  an  adiabatic  wall  boundary  condition  was  used.  Results  for  viscous  first-order 
PPM,  EFMO  and  a  Riemann  solver  due  to  Jacobs  (1991)  are  compared.  The  first-order  EFM 
method  tends  to  thicken  the  boundary  layer  and  underestimates  the  value  of  the  wall-temperature 
as  calculated  with  a  recovery  factor  r  =  1.  There  is  general  good  agreement  between  EFMO  and 
Jacob’s  method  (which  does  not  handle  strong  normal  shocks  weU).  Figure  IV. 6  shows  results  for 
laminar  flow  past  a  cone  at  incidence  24^.  The  other  flow  conditions  are  as  above.  This  case  was 
chosen  to  match  the  experimental  conditions  of  Tracy  (1969).  The  parts  of  Figure  IV. 6  show  the 
Mach-number  distribution  and  the  sectional  streamlines.  The  Mach-number  contours  clearly  show 
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(a) 


(b) 


(c) 


Figure  IV. 6.  Cone  of  half-angle  lO*’  at  incidence  24^  in  a  hypersonic  stream  at  M  =  7.95.  Frozen  flow, 
Re  =  4.2  X  10®  based  on  distance  from  apex,  (a)  Contours  of  Mach  number,  (b)  Sectional 
streamlines  in  a  plane  normal  to  the  cone  axis,  (c)  Detail  of  streamline  pattern  near  the 
leeward  plane  of  symmetry. 

the  separating  shear  layer  and  the  formation  of  a  weak  shock  just  outboard  of  this  layer  on  the 
leeward  side  of  the  body.  Between  the  shear-layer  and  the  body  a  compact  vortex  forms  near  the 
leeward  plane  of  symmetry.  The  angular  position  of  separation  can  be  estimated  from  the  sectional 
stream-line  pattern  and  is  found  to  be  at  27®  from  the  leeward  symmetry  plane  compared  to  24® 
found  by  Tracy  (1969)  from  surface  stream-hne  patterns. 

IV. 3. 3  EFMO  and  the  IDG  Chemistry  model 

Two  versions  of  incorporating  the  Lighthill-Freeman  ideal-dissociating  gas  (IDG)  chem¬ 

istry  model  exist  and  are  currently  undergoing  testing.  This  strategy  has  been  developed  with  the 
aim  of  comparing  the  predictions  of  the  two  codes  against  each  other,  and  against  experiment  and 
other  computation,  for  the  purpose  of  code  validation  The  first  of  these  {PGP3D.3;  Shaun  Shariff) 
is  based  on  an  operator-splitting  approach  in  which  the  hydrodynamic  and  the  chemistry  terms  in 
the  equations  of  motion  are  split,  and  each  is  implemented  separately  at  each  time  step.  The  result 
is  an  exphcit  time-marching  method  for  both  the  gas-dynamics  and  the  chemistry. 

In  the  second  approach,  the  implicit  PGP3D/ EFMO  code  discussed  above  {PGP3D.2)  has  been 
extended  to  include  IDG  chemistry  in  a  partially  implicit  way  {PGP3D.4'-,  Jean-Marc  Moshetta). 
This  is  done,  as  opposed  to  a  fuUy-implict  strongly- coupled  approach,  because  owing  to  the  presence 
of  source  terms  in  the  chemical  rate  equations,  the  Riemann  invariants  in  the  sense  of  Lax  are  no 
longer  defined,  and  it  is  again  necessary  to  decouple  the  gas  dynamics  part  from  the  chemistry 
equation.  In  the  context  of  EFMO  this  then  gives  an  approximate  Riemann  problem  which  must 
be  solved  with  different  values  of  the  specific  heat  ratio  7  for  the  left  state  and  for  the  right  state. 
This  leads  to  a  local  iterative  procedure  in  order  to  obtain  the  intermediate  states.  This  procedure 


Figure  IV. 7.  Two-dimensional  inviscid  flow  of  dissociating  oxygen  past  a  cylinder.  Moo  =  15.0.  Di¬ 
mensionless  reaction-rate  parameter  =  0.218.  (a)  contours  of  density,  (b)  contours  of 
dissociation  fraction  a. 

has  been  modified  to  guarantee  convergence  in  2  or  3  iterations  by  the  method  of  Abgrall  & 
Montagne  (1989).  The  IDG  rate  equation  is  solved  implicitly  by  second-order  Runge-Kutta,  while 
the  gas-dynamics  is  handled  implicitly  as  described  above. 

Both  PGP3D.3  and  PGP3D.4  are  currently  undergoing  testing  for  a  variety  of  one- dimensional 
and  two-dimensional  cases.  These  include  the  simple  one- dimensional  steady  propagating  shock 
wave  and  two-dimensional  chemically  non-equilibrium  flows  about  cylinders  and  spheres.  In  Figure 
IV. 7  we  show  results  from  a  test  calculation  of  the  inviscid  two-dimensional  flow  of  dissociating 
oxygen  at  Moo  =  15  past  a  cylinder  using  PGP3D.3.  The  dimensionless  reaction-rate  parameter 
is  fi  =  0.218  based  on  the  dissociation  rate  da/dt  on  the  downstream  side  of  a  stationary  normal 
shock  at  the  free-stream  conditions.  For  the  viscous  hypervelocity  flow  about  spheres,  detailed 
comparisons  with  the  computations  and  experiments  of  Chihyung  Wen  (Phd,  Caltech  1994)  are 
planned  for  the  near  future. 

IV. 3. 4  Future  Work 

Fully  three-dimensional  tests  of  (PGP3D,2)  on  the  hypersonic  flow  about  the  a  blunt  nosed 
delta  wing  of  70*^  sweep,  at  an  angle  of  attack  of  30^  and  at  Moo  =  8.7  are  currently  underway  with 
the  aim  of  comparing  these  with  the  calculations  of  Mallet  et  al  1994,  done  using  the  viscous  EFM 
method.  For  the  laminar  viscous  flow  two  Reynolds  numbers,  of  Re  =  2.25  X  10^  and  2.25  X  10^ 
will  be  used.  When  the  one-  and  two-dimensional  tests  on  PGP3D.3  and  PGP3D.4  have  been 
completed,  it  is  planned  to  undertake  computations  of  the  blunt  delta-wing  flow  with  active  non- 
equihbrium  gas  chemistry.  It  is  planned  that  Shaun  Shariff  wiU  use  both  PGP3D,3~4  to  perform  a 
sequence  of  computations  on  viscous  reacting  cone  flow  with  a  view  to  investigation  of  the  chemical 
scaling  of  both  the  windward  and  the  leeward  flow  characteristics.  In  the  longer  term  it  is  planned 
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to  incorporate  a  more  complex  partial  equilibrium  model  of  air  (Rein  1992)  into  PGP3D. 
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CHAPTER  V 

DIAGNOSTICS 


Introduction 

Since  the  major  part  of  the  research  in  this  program  is  experimental,  much  of  our  work  is  also 
oriented  toward  development  of  better  diagnostic  methods.  In  some  instances  the  techniques  being 
developed  are  very  closely  coupled  to  particular  projects  and  the  progress  in  their  development  can 
be  more  effectively  presented  together  with  the  progress  in  the  sub-project  in  which  they  are  being 
applied.  This  is  the  case  for  the  development  of  the  new  holographic  interferometer  described  in 
section  II.2,  and  the  image- correlation  velocimetry  technique  described  in  section  III. 

In  the  case  of  the  following  section  on  laser-induced  thermal  acoustics,  the  presentation  is  sep¬ 
arated  from  a  specific  application,  because  it  is  a  project  in  its  own  right.  The  method  is  now 
established  as  a  powerful  diagnostic  tool,  however,  and  is  ready  to  be  applied  to  the  other  sub- 
projects. 
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V.l  DIAGNOSTICS  WITH  LASER-INDUCED  THERMAL  ACOUSTICS 


Objectives  and  Status  of  Research 

The  objectives  of  the  research  described  in  this  subsection  are  to  develop  laser-induced  thermal 
acoustics  (LITA),  a  nonlinear  optical  technique  for  single-shot  measurement  of  gas  and  flow  prop¬ 
erties,  including  species  concentration,  temperature,  transport  properties,  and  velocity,  and  to 
apply  LITA  to  hypersonic  flows  in  T5. 

This  year,  LITA  has  been  transformed  from  a  promising  technique  into  an  accurate  and  well- 
understood  tool  for  the  study  of  gas  properties.  This  transformation  has  been  accomplished  through 
improvements  in  LITA  signal  modeling  and  a  successful  series  of  LITA  experiments.  The  experi¬ 
ments  demonstrate  the  accuracy  of  the  improved  model  and  the  ability  to  extract  accurate  values 
of  the  sound  speed  and  thermal  diffusivity  from  LITA  signals.  The  experiments  also  demonstrate 
the  ability  to  obtain  LITA  signals  in  gases  over  a  wide  pressure  range  (~3  Pa-14  MPa  at  room 
temperature)  and  in  challenging  environments,  e.g.,  CO2  near  its  liquid-vapor  critical  point,  and 
in  flowing,  turbulent  gases. 


V.1.1  Hardware 

A  high-pressure  bomb  was  built  with  optical  access  that  allows  the  study  of  LITA  signals  and  gas 
behavior  at  various  pressures,  temperatures,  and  flow  conditions.  The  bomb  and  its  gas-handling 
system  cover  a  pressure  range  from  ~3  Pa  (absolute)  to  >20  MPa.  The  design  of  the  gas  handling 
system  ehminates  sources  of  contamination  and  exposure  to  the  test  gases.  The  bomb  was  used 
for  all  but  the  initial  LITA  experiments.  Most  of  the  experiments  were  conducted  using  stationary 
gases,  although  the  effects  of  turbulence  and  flow  were  studied  using  an  interior  stirring  fan.  The 
bomb  is  currently  being  outfitted  as  a  test  section  of  a  combustion-driven  shock  tube  (described 
in  Belanger,  1993).  This  shock  tube  will  permit  LITA  studies  of  high-speed  and  high  temperature 
gases  using  the  current  LITA  apparatus,  thereby  easing  the  transition  of  LITA  from  static  ceU  to 
the  T5  free-piston  shock  tunnel. 


V.1.2  Theoretical  model 

Progress  has  also  been  made  on  the  theoretical  modeling  of  LITA  signals.  An  expression  for  the 
LITA  signal  had  previously  been  derived  including  some  of  the  effects  of  the  use  of  spatially  finite 
laser  beams  (Cummings,  1992).  This  expression  successfully  predicted  the  nature  and  amplitude 
of  LITA  signals  from  the  first  round  of  LITA  experiments.  However,  there  were  discrepancies  in 
some  cases  between  the  measured  and  theoretical  signals  that  arose  from  neglected  beam-size  effects 
(Cummings,  1994b).  Therefore,  an  expression  for  the  LITA  signal  that  includes  complete  beam-size 
effects  was  derived.  Through  this  derivation,  phenomenological  interpretations  of  LITA  scattering 
were  elucidated. 

The  analysis  showed  how  two  compelling  but  contradictory  physical  pictures  of  LITA  scattering  were 
correct  in  the  hmits  of  small  and  large  detection  angle  (Cummings  et  al.,  1995a).  The  dependence 
of  LITA  spectra  on  molecular  gas  properties  was  treated  theoretically.  The  effects  of  finite  laser 
noise  and  flow  velocity  were  also  explored.  This  expression  accurately  models  LITA  signals  over 
the  range  of  experimental  conditions,  except  at  the  lowest  pressures,  where  the  theory  needs  to  be 
modified  to  include  rarefied  gas  effects. 
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V.1.3  Comparison  with  experiment 

Examples  of  the  agreement  between  theoretical  and  measured  signals  appear  in  Fig.  V.l  and 
Fig.  V.2.  Fig.  V.l  shows  experimental  and  theoretical  time-histories  of  LITA  signals  formed  by 
the  opto-acoustic  effect  of  thermalization  (Cummings,  1995b)  at  relatively  low  pressure  in  air  at 
~300  K.  seeded  with  trace  (<  10  parts-per-million)  NO2.  The  experimental  apparatus  used  to 
obtain  these  signals  is  discussed  in  Cummings,  1994b  and  Cummings,  1995b.  Signal  averaging 
(3-100  shots)  was  typically  employed  in  these  experiments  so  that  systematic  deviations  between 
experimental  and  theoretical  signals  would  be  unambiguously  exposed.  However,  at  relatively  high 
seed-gas  concentration  or  high  buffer  gas  pressure,  signal  averaging  was  not  strictly  needed  even 
for  the  stringent  task  of  model  verification. 

time  after  driver  pulse  (ns) 
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Fig.  V.l.  Experimental  and  theoretical  LITA  signals  from  thermalization  in  NOo-seeded  air  near- 
atmospheric  pressure  and  300  K. 


The  grating  Knudsen  number — the  ratio  of  the  molecular  mean-free  path  in  the  sample  volume  to 
the  laser-induced  grating  wavelength — is  the  most  important  parameter  in  LITA  scattering.  LITA 
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typically  excels  when  this  parameter  is  much  less  than  unity.  Pig.  V  .2  shows  a  LITA  signal  taken  in 
seeded  air  at  3.45  MPa  and  300  K.  The  long  oscillation-free  tail  of  the  signal  is  caused  by  the  motion 
of  the  laser-induced  acoustic  traveling  waves  away  from  the  source  beam.  A  detailed  discussion  of 
the  physics  from  which  these  signal  patterns  arise  appears  in  Cummings,  1995b,  along  with  plots 
of  the  experimental  and  theoretical  signals  over  a  wider  range  of  grating  Knudsen  numbers. 

time  after  driver  pulse  (ns) 

-1000  1000  3000  5000  7000  9000 


Fig.  V.2.  E.xperimental  and  theoretical  LITA  signal  from  thermalization  in  N02-seeded  air  at  3.45 
MPa  and  300  K. 


In  gases  with  low  levels  of  seeding  or  strong  nonresonant  susceptibility,  the  opto-acoustic  effect 
of  electrostriction  (Cummings  et  al..  1995c)  contributes  to  the  LITA  signal.  .A  combination  of 
electrostriction  and  thermalization  produced  the  signal  that  appears  with  a  fit  theoretical  signal  in 
Fig.  V.3,  taken  in  unseeded  atmospheric-pressure  CO2  at  .300  K. 

Automated  fitting  of  theoretical  and  experimental  curves  was  performed  using  a  nonlinear  least- 
squares  numerical  technique.  This  fitting  is  a  type  of  optimal  nonlinear  filtering  which  permits 
accurate  measurements  of  signal  parameters,  e.g..  sound  speed  and  thermal  diffusivity.  from  noisy 
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time  after  driver  pulse  (ns) 


Fig.  V.3.  Experimental  and  theoretical  LITA  signals  from  electrostriction  and  thermalization  in 
unseeded  CO2  at  atmospheric  pressure  and  300  K.  The  secondary  peaks  in  this  signal 
arise  from  electrostriction. 

signals.  Currently,  the  automated  filtering  of  a  2000-point  LITA-signal  time  history  requires  typi¬ 
cally  one  minute  on  a  Sun  Sparc  5.  An  example  of  sound-speeds  extracted  from  LITA  signals  taken 
in  N02-seeded  air  at  different  pressures  and  ~300  K  appears  in  Fig.  V.4. 

The  sound  speeds  are  corrected  for  small  temperature  excursions  by  scaling  them  with  i/300K/r, 
where  T  is  the  measured  temperature  of  the  gas  in  the  test  ceU,  and  normalizing  by  their  proper 
value  at  atmospheric  pressure.  The  data  of  Hilsenrath  et  al.,  (1955),  appear  as  solid  circles.  The 
agreement  between  the  measured  and  published  values  is  ~0.1%  when  the  pressure  is  above  ~200 
kPa.  At  lower  pressures.  finite-Knudsen-number  effects  and  rotational  nonequilibrium  effects  arise, 
hmiting  the  accuracy  of  the  current  model,  which  is  based  upon  the  equations  of  hydrodynamics. 

Fig.  V.5,  a  plot  of  the  thermal  conductivity  (derived  from  the  thermal  diffusivity  extracted  from 
LIT.A.  signals)  vs.  pressure,  demonstrates  several  shortcomings  of  the  current  theoretical  model. 
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Fig.  V.4.  Sound  speed  of  air  vs.  pressure  corrected  for  300  K  and  normalized  to  the  value  at 
atmospheric  pressure.  The  solid  circles  are  independent  measurements  from  the  National 
Bureau  of  Standards. 

Two  data  sets  taken  with  slightly  different  beam  geometries  are  compared  with  the  data  of  Var- 
gaftik.  et  ah,  (1994).  At  low  pressures,  the  model’s  treatment  of  the  thermalization  process  is 
oversimplified,  ignoring  long-time-scale  deposition  of  excited-state  energy  into  thermal  energy  of 
the  gas.  In  experiments  this  slow  deposition  prolongs  the  life  of  the  thermal  grating.  The  auto¬ 
mated  filter  reacts  by  artificially  lowering  the  extracted  thermal  diffusivity.  An  improved  treatment 
of  thermalization  dynamics  in  the  LITA  signal  model  can  practically  eliminate  this  effect,  although 
such  a  treatment  could  require  extensive  study  of  the  quenching  behavior  of  the  seed  species  in  the 
buffer  gas. 

At  moderate  pressures,  the  effects  of  slight  beam  misalignment  skew  the  shape  of  the  signal  in  a 
manner  that  is  consistent  with  an  decreased  thermal  diffusivity.  In  one  data  series,  the  data  rise 
rapidly  at  high  pressure  to  within  ~1%  of  the  published  value  after  fine  adjustments  were  made  in 
the  beams  in  order  to  optimize  their  overlap.  As  the  pressure  was  reduced,  differential  refraction 
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Fig.  V.5.  Thermal  conductivity  of  air  from  LITA  measurements  and  data  of  Vargaftik  et  ai. 

again  caused  the  beam  alignment  to  degrade.  While  beam  misalignment  effects  could  be  trivially 
inserted  into  the  current  model,  a  better  solution  to  this  modeling  error  in  future  experiments  could 
be  to  employ  unfocused  beams,  which  would  dramatically  ease  alignment  requirements. 

The  LITA  model  predicts  and  experiments  verify  that  the  peak  LITA  signal  intensity  scales  as 
the  square  of  the  buffer-gas  pressure  at  constant  seed-gas  concentration,  provided  the  seed-gas 
absorption  spectrum  is  relatively  constant  with  pressure.  The  integrated  signal  intensity  scales  as 
the  cube  of  the  buffer-gas  pressure.  Fig.  V.6  shows  the  pressure  dependence  of  peak  and  integrated 
signal  intensities  from  experiments.  The  approximate  regimes  over  which  various  complications 
may  arise  in  taking  and  interpreting  LITA  signals  are  indicated. 

The  effects  of  convection  and  turbulence  in  the  sample  volume  on  LITA  signals  were  observed 
and  analyzed.  The  magnitude  of  these  effects  now  may  be  predicted  from  general  flow  properties, 
such  as  the  Mach  and  Reynolds  numbers  and  beam  parameters,  such  as  the  frequencies,  crossing 
angles,  and  diameters  (Cummings,  1995b).  Flow-related  difficulties  may  be  mitigated  in  advance 
of  experiments  by  adjusting  experimental  parameters  such  as  the  beam  geometry. 

V.1.4  Plans 

A  new  series  of  LITA  tests  will  employ  the  combustion-driven  shock  tube  throughout  the  summer 
and  into  the  fall.  In  these  and  a  continued  series  of  static  tests,  new  approaches  to  LITA  will  be 
explored,  including  the  use  of 

•  advanced  techniques  to  ease  beam-alignment. 


•  ultra-violet  light  to  study  LITA  using  the  target  species  NO,  a  gas  indigenous  to  T5, 

•  a  flashlamp- pumped  dye  laser  as  a  source  laser  to  boost  signal  levels, 

•  heterodyne  detection  to  recover  the  Doppler  shift  of  the  signal  beam  scattered  by  convecting 
gratings,  thus  measuring  the  gas  velocity 

and  a  variety  of  other  approaches  designed  to  reduce  the  risk  and  increase  the  rewards  of  LITA 
measurements  in  T5. 
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Interactions 

The  success  of  this  project  has  led  to  intensive  interest  from  the  company  MetroLaser.  Two 
contracts  with  MetroLaser  have  resulted.  Under  these  contracts  MetroLaser  bring  equipment  and 
instrumentation  to  the  T5  laboratory,  and  use  T5-lab  equipment  and  instrumentation  to  further 
develop  LITA,  and  to  apply  resonant  holographic  interferometric  spectroscopy  (RHIS)  in  the  T5 
environment. 
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Presentations  at  GALCIT  Fluid  Mechanics  Seminar,  and  at  GALCIT  Research  Conference.  Sem¬ 
inar  presentation  at  Princeton  University.  Discussions  with  the  group  of  Stepanek  at  AEDC. 
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